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1 Introduction 



The mechanism for confinement of quarks and gluons into colorless hadrons still is 
little understood. Some theoretical insight could be obtained from disproving the 
cluster decomposition property for color-nonsinglet gauge-covariant operators. One 
idea in this direction is based on the possible existence of severe infrared divergences, 
i.e., divergences which cannot be removed from physical cross sections by a suitable 
summation over degenerate states by virtue of the Kinoshita-Lee-Nauenberg theorem 
Q. Such divergences can provide damping factors for the emission of colored states 
from color-singlet states (see Q). However, the Kinoshita-Lee-Nauenberg theorem 
applies to non-Abelian gauge theories in four dimensions order by order in perturba- 
tion theory 1^] . Therefore, such a description of confinement in terms of perturbation 
theory is impossible. In fact, extended to Green's functions, the absence of unphysical 
infrared divergences implies that the spectrum of QCD necessarily includes colored 
quark and gluon states to every order in perturbation theory [Q]. 

An alternative way to understand the insufficiency of perturbation theory to ac- 
count for confinement in four dimensional field theories is that confinement requires 
the dynamical generation of a physical mass scale. In presence of such a mass scale, 
however, the renormalization group [RG] equations imply the existence of essential 
singularities in physical quantities, such as the 5-matrix, as functions of the cou- 
pling at 5 = 0. This is because the dependence of the RG invariant confinement scale 
on the coupling and the renormalization scale ^ near the ultraviolet fixed point is 
determined by [|| 

A = ^exp(-/'^)^^Vexp(-^), /3o > 0. (1) 

Since all RG invariant masses in massless QCD will exhibit the behavior (|^) up to a 
multiplicative constant the ratios of all bound state masses are, at least in the chiral 
limit, determined independent of all parameters. 

Therefore, to study the infrared behavior of QCD amplitudes non-perturbative 
methods are required. In addition, as singularities are anticipated, a formulation in 
the continuum is desirable. One promising approach to non-perturbative phenomena 
in QCD is provided by studies of truncated systems of its Dyson-Schwinger equa- 
tions [DSE], the equations of motion for QCD Green's functions. Typical truncation 
schemes resort to additional sources of information like the Slavnov-Taylor identities, 
as entailed by gauge invariance, to express vertex functions and higher n-point func- 
tions in terms of the elementary two-point functions, i.e., the quark, ghost and gluon 
propagators. In principle, these propagators can then be obtained as selfconsistent 
solutions to the non-linear integral equations representing the closed set of truncated 
DSEs. 

The underlying conjecture to justify such a truncation of the originally infinite 
set of DSEs is that a successive inclusion of higher n-point functions in selfconsistent 
calculations will not result in dramatic changes to previously obtained lower n-point 
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functions. To achieve this it is important to incorporate as much independent infor- 
mation as possible in constructing those n-point functions which close the system. 
Such information, e.g., from implications of gauge invariance or symmetry properties, 
can be relied on being reproduced from solutions to subsequent truncation schemes. 

So far, available solutions to truncated Dyson-Schwinger equations of QCD do 
not even fully include all contributions of the propagators itself. In particular, even 
in absence of quarks, solutions for the gluon propagator in Landau gauge used to 
rely on neglecting ghost contributions ^, ^, ^ which, though numerically small in 
perturbation theory, are unavoidable in this gauge. While this particular problem 
can be avoided by ghost free gauges such as the axial gauge, in studies of the gluon 



Dyson-Schwinger equation in the axial gauge |1C, |ll], 12], the possible occurrence of 
an independent second term in the tensor structure of the gluon propagator has been 
disregarded [|^]. In fact, if the complete tensor structure of the gluon propagator 
in axial gauge is taken into account properly, one arrives at a coupled system of 
equations which is of similar complexity as the ghost-gluon system in the Landau 
gauge and which is yet to be solved |0 . 

In addition to providing a better understanding of confinement based on studies 
of the behavior of QCD Green's functions in the infrared, DSEs have proven suc- 
cessful in developing a hadron phenomenology which interpolates smoothly between 
the infrared (non-perturbative) and the ultraviolet (perturbative) regime, for recent 
reviews see, e.g., IC]. In particular, a dynamical description of the spontaneous 
breaking of chiral symmetry from studies of the DSE for the quark propagator is 
well established in a variety of models for the gluonic interactions of quarks [17|. For 
a sufficiently large low-energy quark-quark interaction quark masses are generated 
dynamically in the quark DSE in some analogy to the gap equation in superconduc- 
tivity. This in turn leads naturally to the Goldstone nature of the pion and explains 
the smallness of its mass as compared to all other hadrons |18|. In this framework 
a description of the different types of mesons is obtained from Bethe-Salpeter equa- 
tions for quark-antiquark bound states |19|. Recent progress towards a solution of a 
fully relativistic three-body equation extends this consistent framework to baryonic 



bound states |g, P, |2|. 



In this paper we present a truncation scheme for the Dyson-Schwinger equations 
of QCD in Landau gauge which implements the Slavnov-Taylor identities for the 
3-point vertex functions. Neglecting contributions from 4-point correlations such as 
the 4-gluon vertex function and irreducible scattering kernels, this yields a closed set 
of non-linear integral equations for the propagators of gluons, ghosts and, in an obvi- 
ous extension, also quarks. We present a solution obtained for the pure gauge theory 
without quarks in a combination of numerical and analytical methods. The organi- 
zation of the paper is as follows: The general framework is outlined in section 2. In 
section 3 we review some earlier results in a scheme proposed by Mandelstam which 
is most easily introduced by further simplifying the general scheme of the preceeding 
section, and we comment on available axial gauge studies of the gluon propagator. 
We concentrate on the couped system of DSEs for gluons and ghosts from section 4 



2 



on, in which we introduce an approximation to reduce this system to a set of one- 
dimensional non-hnear integral equations. This approximation is especially designed 
to preserve the leading behavior of these propagators in the infrared, and we give an 
analytic discussion of this behavior. In section 5 we introduce the non-perturbative 
momentum subtraction scheme and discuss the accordingly extended definition of 
the corresponding running coupling. The actual details of the renormalization, the 
definition of the renormalization constants and the renormalized ultraviolet finite 
equations are given in section 6. An extended analytic discussion of the solutions to 
the renormalized equations in terms of asymptotic expansions, a necessary prerequi- 
site to a numerical solution, is presented in section 7. The numerical results and their 
discussion will be the issue of section 8. The non-perturbative result for the running 
coupling, i.e., the strong coupling as function of the renormalization scale, and the 
corresponding Callan-Symanzik /3-function are given in section 9. We compare our 
results to recent lattice calculations of the gluon and ghost propagators in section 10, 
discuss the implications on confinement for gluons in section 11, and we wrap things 
up with our summary and conclusions in section 12. We wish to point out that one 
main result, the existence of an infrared stable fixed point at a critical value of the 
coupling ac — 9.5, which could, in principle, be obtained independent of the details 
of the renormalization as well as the numerical methods and results, is found in our 
previous publication |23]. This result, meaningless in absence of a complete solution, 
will be incorporated also here in order to make the present paper self-contained. 



2 A Truncation Scheme for the Propagators of QCD in 
Landau Gauge 

Besides all elementary two-point functions, i.e., the quark, ghost and gluon propaga- 
tors, the Dyson-Schwinger equation for the gluon propagator also involves the three- 
and four-point vertex functions which obey their own Dyson-Schwinger equations. 
These equations involve successive higher n-point functions, see figure ||. For simplic- 
ity we consider the pure gauge theory and neglect all quark contributions. As for the 
truncation scheme, the extensions necessary to include quarks in future calculations 
are straightforward and will be given at the end of this section. 

A first step towards a truncation of the gluon equation is to neglect all terms 
with four-gluon vertices. These are the momentum independent tadpole term, an 
irrelevant constant which vanishes perturbatively in Landau gauge, and explicit two- 
loop contributions to the gluon DSE. The latter are subdominant in the ultraviolet, 
and we can therefore expect solutions to have the correct behavior for asymptotically 
high momenta without those terms. In the infrared it has been argued that the 
singularity structure of the two-loop terms does not interfere with the one-loop 
terms [^^. It therefore seems reasonable to study the non-perturbative behavior of 
the gluon propagator in the infrared without explicit contributions from four-gluon 
vertices, i.e., the two-loop diagrams in figure |^. Without quarks the renormalized 
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Figure 1: Diagrammatic representation of the gluon, gliost and quark Dyson- 
Schwinger equations of QCD. The wiggly, dashed and sohd hnes represent the prop- 
agation of gluons, ghosts and quarks respectively. A filled blob represents a full 
propagator and a circle indicates a one-particle irreducible vertex. 
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Figure 2: Diagrammatic representation of the gluon and ghost Dyson-Schwinger 
equations of QCD without quarks. In the gluon DSE terms with four-gluon vertices 
have been dismissed. 

equation for the inverse gluon propagator in Euclidean momentum space with positive 
definite metric, g^^y = (5^i,, (color indices suppressed) is thus given by 

D-lik) = Z^D%l{k) - g^N.Zr J ^ iq^ Dg{p) Doiq) GAq,p) (2) 

+ g'^NcZi-J ■j^--^r^p^{k,-p,q)Dap{q)Dp^{p)T,3„^{-q,p,-k) , 

where p = k + q, D*' and F" are the tree-level propagator and three-gluon vertex, Dq 
is the ghost propagator and F and G are the fully dressed 3-point vertex functions. 
The DSE for the ghost propagator in Landau gauge QCD, without any truncations, 
will remain unchanged in its form when quarks are included, 

Dc\k) = -Zs k" + g^N, Zi J ik^ D^,{k - q) G,{k, q) Dciq) ■ (3) 

The coupled set of equations for the gluon and ghost propagator, eqs. @ and (^), is 
graphically depicted in figure |^. The renormalized propagators for ghosts and gluons, 
Dg and Dfj_y, and the renormalized coupling g are defined from the respective bare 
quantities by introducing multiplicative renormalization constants, 

Z3L»g:=^G> Z^D^,:=D%, Zgg := go . (4) 

3/2 1/2 — 

Furthermore, Zi = ZgZ-^', Zi = ZgZ^'^Zs. We use the following notations to sepa- 
rate the structure constants f"''"^ of the gauge group SU {Nc = 3) (and the coupling 
g) from the 3-point vertex functions: 

M-, a ^ 

TX{k,P,q)=9r''{27r)^SHk + P + q)T^,pik,p,q) . A • (5) 
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The arguments of the 3-gluon vertex denote the three incoming gluon momenta 
according to its Lorentz indices (counter clockwise starting from the dot). With this 
definition, the tree-level vertex has the form, 

r)!i.p(^,P,9) = -i{k-p)p5f,y - i{p - q)^5yp - i{q - k)^5pp . (6) 

The arguments of the ghost-gluon vertex are the outgoing and incoming ghost mo- 
menta respectively, 

k- 

Gf%p, q) = gr'^Gpip, q) = gr'Hq^G^AP, q), a 

X p 



where the tensor G^i,{p,q) contains the ghost-gluon scattering kernel (for its defi- 
nition see, e.g., ref. ||2^). At tree-level this kernel vanishes which corresponds to 
GfjLuiPiq) = ^fiu- Note that the color structure of all three loop diagrams in fig. |2| 
is simply given by f^-^^.f^'^'^ = —jsf^^^'b -^j^jgj^ -^g^g in eqs. (|2|,^) suppressing the 
trivial color structure of the propagators ~ 

The ghost and gluon propagators in the covariant gauge, introducing the gauge 
parameter ^, are parameterized by their respective renormalization functions G and 
Z, 

Dcik) = and (8) 

D,Ak) = + 



In order to arrive at a closed set of equations for the functions G and Z, we still have 
to specify a form for the ghost-gluon and the 3-gluon vertex functions. We will now 
propose a construction of these vertex functions based on Slavnov-Taylor identities 
as entailed by BRS invariance. For the 3-gluon vertex the general procedure in 



such a construction was outlined in the literature ||2^, 26, |2^. Since this procedure 
involves unknown contributions from the ghost-gluon scattering kernel, it cannot 
be straightforwardly applied to serve the present purpose, which is to express the 
vertex functions also in terms of the functions G and Z. To this end we have to make 
additional assumptions on the scattering kernel. Since it is related to the ghost-gluon 
vertex, we start with the construction of the latter. 

For our actual solutions later on, we will be using the Landau gauge (^ = 0) in 



which one has Zi = 1 |28|. Therefore, in the case of the ghost-gluon vertex the choice 
of its tree-level form in the DSEs may seem justified to the extend that it would at 
least not affect the leading asymptotic behavior of the solutions in the ultraviolet. 
This is in contrast to the other vertex functions in QCD which have to be dressed at 
least in such a way as to account for their anomalous dimensions, if the solutions to 
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the DSEs in terms of propagators are expected to resemble their leading perturbative 
behavior at short distances, i.e. their correct anomalous dimensions. 

A good way to achieve this is to construct vertex functions from their respective 
Slavnov-Taylor identities. Those identities generally fail, however, to fully constrain 
the vertex functions. Methods of refinement of such constructions from implications 
of multiplicative renormalization have been developed for the fermion-photon vertex 
in QED Our particular focus in the present context is on the 3-gluon vertex 

which has a considerably higher symmetry than fermion vertices. This symmetry 
alleviates the problem of unconstraint terms and we do not expect our results to be 
overly sensitive to such terms in the 3~gluon vertex. 

The Slavnov-Taylor identity for the 3-gluon vertex in momentum space is given 

by 



ikpTpyp{p,q,k) = G{k ) yGp„{-k,q)Vau{q) ^j-^ (9) 

2 ^ 

- G^^{-k,p)Vap{p) ^ 



Z{p^) j ' 

where Vp,u{k) = 5py — kpky/k'^ is the transversal projector. A simple solution to (P) 
is possible if ghosts are neglected completely, i.e., for G{k'^) = 1 and Gpy = 5py (see 
the next section). This is obviously not satisfactory for our present study addressing 
ghost contributions to the gluon propagator in particular. Given that the anomalous 
dimension of the ghost-gluon vertex vanishes in Landau gauge, it may seem appealing 
to retain its tree-level form by setting Gp^ = 5pu while taking the presence of the 
ghost renormalization function G{k'^) into account in solving (|9|). Unfortunately, this 
leads to a contradiction. The Slavnov-Taylor identity for the 3-gluon vertex can be 
shown to have no solution with this Ansatz. One possibility to resolve this problem 
is to add a term corresponding to a non-vanishing scattering kernel of the form, 

Gpu{p,q) = V +o,{p^^q^\k^){^tiuvq - p^qti) ■ (lo) 

While the additional term in (|lO|) corresponding to a non-trivial ghost-gluon scat- 
tering kernel does not contribute to the ghost-gluon vertex (being transverse in g^,), 
it is possible to make a minimal Ansatz for the unknown function a(x, y; z) such 
that (^) can be solved. The proof that (^) cannot be solved for Gpu = 5py and the 
construction of a(x, y; z) are given in appendix ^ This procedure is certainly not 
unambiguous, and it is furthermore against our original intention to neglect explicit 
contributions from irreducible four-point functions as it implies a very specific form 
for the ghost-gluon scattering kernel. A truncation in which the ghost-gluon vertex 
function is replaced with its tree-level form, Gp{p,q) — > iqp, while compatible with 
the desired short distance behavior of the solutions, implies that some non-trivial as- 
sumptions on the otherwise unknown ghost-gluon scattering kernel have to be made 
in order to solve the Slavnov-Taylor identity for the 3-gluon vertex. 
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In light of this, we prefer a second possibihty to resolve the present problem, which 
is to reconsider the form of the ghost-gluon vertex to start the truncation scheme 
with. A suitable source of information is provided by a Slavnov-Taylor identity which 
involves the ghost-gluon vertex directly, similar to the one for the quark-gluon vertex. 
To our knowledge, such an identity has not been stated in the literature so far. We 
will therefore outline its derivation from the usual BRS invariance briefly in the 
following. Neglecting irreducible ghost-ghost scattering, we will arrive at an identity, 
which allows to express the ghost-gluon vertex in terms of the ghost renormalization 
function G(/c^). The vertex constructed this way will have the correct short distance 
behavior (no anomalous dimension) and it will constrain G^y{p,q) just enough to 
admit a simple solution to the 3-gluon Slavnov-Taylor identity @ without further 
assumptions. 

In order to construct a non-perturbatively dressed ghost-gluon vertex, we start 
from BRS invariance of the generating functional of pure Yang-Mills theory, 

5^^sZ%v,J] = ^ (11) 
A[??, r?, J] := l^j d^x (^rD'^'^c^ - ^p^^f)''c^c^ - ^dA'^r^''^ ^ _ = , 

where D""^ = 5°'^d+gf°'^'^A'^ is the covariant derivative. We retained non-zero sources 
{??, r/, J} for the ghost and gauge fields {c, c. A} in order to perform suitable deriva- 
tives. Fermion fields have been omitted for simplicity. Their inclusion would not 
affect the present derivation. We find that the particular derivative of (11) with 



respect to the sources, suitable to constrain the ghost-gluon vertex, is given by: 
6^ 



5rj^{z)5rf'{y)5ri"-{x] 



A[r],r],J] 



= , (12) 

ri=ri=J=0 



and we obtain 

Uc'{z)c^{y)dA<'ix)) - Uc''{z)c^{x)dA\y)) (13) 

= -p'"''{c''{z)c%z)(^{x)c^{y)) . 

Note that this symbolic notation refers to full reducible correlation functions. The two 
terms on the l.h.s. of ( p^ ) can be decomposed in the ghost-gluon proper vertex and 
the respective propagators. The derivative on the gluon leg thereby projects out the 
longitudinal part of the gluon propagator which, by virtue of its own Slavnov-Taylor 
identity k^D^i,{k) = ^ k^/k"^, remains undressed in the covariant gauge (see eq. (^)). 
The r.h.s. contains a disconnected part plus terms due to ghost-ghost scattering. 
At this point we make the truncating assumption mentioned previously. We neglect 
ghost-ghost scattering contributions to the vertex and retain only the reducible part 
of the correlation function on the r.h.s. of (|l^) corresponding to disconnected ghost 
propagation, i.e., 

(c^c^cV) = {c''c^){c''c^) - {c''^) {c'c^) + connected. (14) 
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This assumption results in an Abelian Ward-Takahashi type of equation for the 
ghost-gluon vertex, relating its longitudinal part to a sum of inverse ghost propaga- 
tors. In momentum space with the ghost renormalization function as introduced in 
(P) we obtain (for details see appendix |A|), 

ik,Gf%p, q) - Gik^) ippGf%-k, q) = gP^' ^'^g(g2^^^'^ • (15) 

Note that this equation is valid for all ^ including ^ = (Landau gauge) . Eq. (|l^) is 
solved by the following particularly simple form for the vertex which was also used 



in ref. |3g|l, 



Gf''{p,q)=gr'"iq^.^^ ■ (16) 



G{q^ 

It is important to note, however, that the solution ( [T^ ) to eq. ( [T5| ) violates the 
symmetry properties of the vertex explicitly. In order to check, whether a symmetric 
solution to (p!5|) exists, we make the general Ansatz, 



G^{p,q) = i{p + q)^,a{k^;p^,q^) - i{p - q)^,h{k^]P^ ,q^) • (17) 

The reason for this particular choice is that in general covariant gauges one has 

a{x]y,z) = a{x]z,y) , (18) 

whereas the symmetry of h{x\ y, z) is undetermined. In Landau gauge, however, the 
mixed symmetry term proportional to b does not contribute to observables directly, 
since it is purely longitudinal in the gluon momentum. With the Ansatz (|l7|) we find: 



G{x) ^(z - x) a{y; x,z)+y b{y; x, z)^ (19) 
+ G{y)(^{z-y)a{x;y,z) + xb{x;y,z)^ = ^^(^^'^(y) 



G{z) 



From all terms proportional to z it follows that 



a{x;y,z) = ^ ^rr^ + ai{x,y; z)-^ , ai{x,y;z) = -ai{y,x;z) . 
2 G{z) G{y) 

Here the undetermined antisymmetric function oi(x,y;2;) has to guarantee the sym- 
metry of a{x;y,z) c.f. ([l8|) . This can be achieved most easily by setting 

1 fGix) ^ G{x) 



ai{x,y]z) = ^ {g{x) - G{y)^ a{x;y,z) 



2\G{z) G{y) 

The mixed symmetry term follows readily from terms proportional to x or y in (|l^). 



N G(x) , , 1 /G'(x) G{x) , 
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Collecting the above terms, from (IT^ the ghost-gluon vertex can be written, 



= + ^py - 1 j + longitudmal terms . (22) 

In the last line we rearranged the terms and displayed explicitly only the ones relevant 
for the vertex in Landau gauge. While terms longitudinal in the gluon momentum 
k = p — q are irrelevant in the vertex, they contribute to the 3-gluon vertex Slavnov- 
Taylor identit y (^ ). The form of G^MuiPiQ) necessary to generate the complete ghost- 
gluon vertex (|2l|), is given by 



G{k^) ^(G{k^) \p,q. 



Here, the last term does not contribute to (P). Therefore, while the relevant term 
to the ghost-gluon vertex in Landau gauge is given by the one explicitly stated in 
(p^), the term of the tensor (^) to be used in dD is G^,{p,q) = (G(A;2)/G((72)) V 
and not what would follow from ( [2^ ) with the longitudinal terms neglected before 
constructing G, i.e., {G{k'^)/G{q^)+G{k'^)/G{p'^)-l)5^y. Additional contributions to 
G/iuip, q) which are purely transverse in q^ can arise from the ghost-gluon scattering 
kernel. Such terms cannot be constraint from the form of the ghost-gluon vertex. In 
contrast to the scheme with the tree-level ghost-gluon vertex, where precisely these 
terms are necessary to solve the 3-gluon Slavnov-Taylor identity, using the dressed 
vertex (^Tj) and (|23|) , unknown contributions from the ghost-gluon scattering kernel 
can be neglected in (^. Thus using ( p3D in (^) we obtain 



ik,T,^,{p, q, k) = Gik') { -PM-^^^^ - 'P^'^iP) ct^lzll^ 1 ' (2^) 



which is only slightly generalized as compared to the Abelian like identity obtained 
neglecting ghosts completely. Using the symmetry of the vertex the solution to ( p^ ) 
fixes the vertex up to completely transverse parts. It can be derived straightforwardly 



along the lines of the general procedure outlined in 27]. In the present case, the 
solution is 

^t^up{p,q,k) = -A+{p^,q^;k^)5^ui{p- q)p - A^{p'^,q^;k^) 6^,yi{p + q)p 

~ 2 ^2 ' q2 ^^i^^PI ~ Piy1tJ.)'^iP ~ Q)p + cyclic permutations , (25) 

with 



AJ,W-.k^) = Cie)U^£^^^^^^] . (26) 
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At this point, additional unconstrained transverse terms in the vertex are ignored. 
A counter example for this to be justified is the fermion-photon vertex in QED, 
for which it has been shown that the transverse part is crucial for multiplicative 
renormalizability (see, refs. ||2^). A similar construction to fix transverse pieces of 
the vertices in QCD is still lacking. The most important difference, however, which 
tends to alleviate the problem of undetermined transverse terms in the case of the 
3-gluon vertex as constructed from its Slavnov- Taylor identity is its high symmetry. 
In contrast to the QED fermion vertex, in which 8 of the 12 independent tensor 
terms are transverse in the photon momentum, in the case of the 3-gluon vertex only 
4 out of 14 terms (involving 2 unknown functions) are transverse in all three gluon 
momenta and thus unconstraint by the Slavnov-Taylor identity (see [^). Therefore, 
its Slavnov-Taylor identity together with its full Bose (exchange) symmetry puts 
much tighter constraints on the 3-gluon vertex than the Ward-Takahashi/Slavnov- 
Taylor identities do on the fermion-photon/gluon vertices. Note also that there are 
no undetermined transverse terms in the ghost-gluon vertex of the present truncation 
scheme. 

A further important difference between fermion vertices and the 3-gluon vertex 
function is that transverse terms in the vertices of (electrons) quarks are of particular 
importance due to their coupling to transverse (photons) gluons in the Landau gauge. 
In contrast, for the present purpose of the 3-gluon vertex function, to be used in 
truncated gluon DSEs, it is well known that the longitudinal contribution to the 3- 



gluon loop is the one of particular importance |31|. This implies that even though the 
two gluons within the loop are transverse in Landau gauge, the third (external) leg 
of the 3-gluon vertex must not be connected to a transverse projector (see below). 
This makes the important difference, because only the terms that are transverse with 
respect to all three gluon momenta are unconstraint in the vertex. 

Now, we have set up a closed system of equations for the renormalization functions 
G{k'^) and Z{k'^) of ghosts and gluons, consisting of their respective DSEs (|2|) and (^) 



with the vertex functions given by ( |22D and (25/^). Thereby we neglected explicit 4 



gluon vertices (in the gluon DSE (0)), irreducible 4-ghost correlations (in the identity 



for the ghost-gluon vertex (15)) and contributions from the ghost-gluon scattering 
kernel (to the Slavnov-Taylor identity (^)). This is the most important guideline 
for our truncation scheme according to the general idea of successively taking higher 
n-point functions into account. 

It is straightforward to extend the present scheme to selfconsistently include the 
quark propagator which has the general structure, 

Sip) = 2^ — ■ (27) 

For the quark-gluon vertex r^(p, q) with momentum arguments p and q for outgoing 
and incoming quarks respectively, the Slavnov-Taylor identity is [32|, 



G-\k^)ik^,Tl{p,q) = (gt'^-B'^{k,q))iS-\p) (28) 
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iS-'{q)lgt''-B''{k,q)] , k=p-q. 



Here, t°' is the SU{3) generator in the fundamental representation, and B°'{k, q) is the 
ghost-quark scattering kernel which again represents irreducible 4-point correlations. 
In the present truncation scheme in which such 4-point correlations are consistently 
neglected, the solution to this Slavnov-Taylor identity is obtained from a particularly 
simple extension to the construction of Ball and Chiu for the solution to the analogous 



Ward-Takahashi identity of Abelian gauge theory [33|. For B"-(k, q) = the Slavnov- 



Taylor identity for the quark-gluon vertex ( pq ) is solved by 

r;ip,q) = -9t'^G{k')\l{A{p')+A{q'))n, + ^±^({A{p') 



^(g^)) o - (b{p^) - B{q^)) ] ) + transverse terms . (29) 




This is justified for a study with emphasis on the infrared behavior of the propa- 
gators even beyond the present truncation scheme, because it has been shown that 
B"'{k,q) — > for A: — > in Landau gauge Note that the only difference to the 
Abelian case considered in at this point is the presence of the ghost renormal- 
ization function appearing in the quark-gluon vertex in Landau gauge. This being a 
minor modification to the structure of the vertex, the ghost renormalization function 
in (^) is, however, crucial for the vertex to resemble its correct anomalous dimension 
in the perturbative limit. Furthermore, our present solution demonstrates that the 
ghost renormalization function G{k'^), being infrared enhanced as we will see below, 
can give an essential contribution to the effective interaction of quarks as part of the 
quark-gluon vertex also in the infrared. 

With the form (|2^ ) for the vertex function in the quark DSE, possibly improved 
by additional transverse terms as discussed for quenched QED in |^, the present 
truncation scheme is extended to a closed set of equations for all propagators of 
Landau gauge QCD as parameterized by the four functions Z{k'^),G{k'^), A(k'^) and 
B{k'^). Its solution would represent for the first time a systematic and complete 
solution to the DSEs of QCD at the level of propagators. In the following, we will 
present our solution to the subset of equations restricted to the pure gauge theory 
without quarks. The selfconsistent inclusion of the quark DSE will be subject to 
further investigations. 



3 The Infrared Behavior of the Gluon Propagator in 
Previous DSE Studies 

Before we turn to the discussion of the coupled system of DSEs for gluons and ghosts 
and its simultaneous solution in an one-dimensional approximation in the following 
sections, we briefly review the current status in DSE studies of the gluon propagator 
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in Landau gauge which is of particular relevance to our present scheme going beyond 
an approximation by Mandelstam that was previously employed in these studies. For 
comparison we also summarize available axial gauge studies of the gluon propagator. 
Their review is necessarily incomplete and we refer the reader to the quoted literature 
for further information. 

A first approximation scheme for the gluon Dyson-Schwinger equation in Lan- 
dau gauge was originally proposed by Mandelstam Q. Compared to the scheme 
outlined above, it consists of a further truncating assumption which is, even though 
working in Landau gauge, to neglect all ghost contributions to the gluon DSE in 
pure QCD (without quarks). As a justification for this, it was usually referred to 
perturbative calculations which yield numerically small ghost contributions to the 
gluon self-energy. Even though there was never any doubt about the importance of 
ghosts for fundamental reasons such as transversality of the gluon propagator and 
unitarity, it was asserted that their quantitive contributions to many hadronic ob- 
servables might remain negligible even beyond perturbation theory. We will see in 
the following sections that our present solutions to the coupled system of gluon and 
ghost DSEs, in the systematic truncation scheme outlined above, yield qualitatively 
quite different results as compared to the Mandelstam approximation, and are thus 
counter examples to this assertion. 

Without ghosts, the solution to the Slavnov-Taylor identity for the 3-gluon vertex 



is obtained from eq. (p3) by setting G = 1, i.e., eq. (25) with 



A,(/,,^t^) = ^,(/,,=) = i(^±^) . (30) 

Assuming that ZijP') is a slowly varying function one may then approximate the 
corresponding solution (p5| ) by 

T^,p{p,q,k) = A+{p',q^)Tl,p{p,q,k) . (31) 

While this form for the full three-gluon vertex simplifies the 3-gluon loop in the gluon 
DSE even more than the use of another bare vertex, i.e., F = F", it was observed by 
Mandelstam to be superior to the latter since it accounts for some of the dressing of 
the vertex as it results from the corresponding Slavnov-Taylor identity. The nature of 
this dressing is such that it cancels the dressing of one of the gluon propagators in the 
3-gluon loop, and without ghost contributions the gluon DSE (^) in the Mandelstam 
approximation thus simplifies to 

D-^ik) = Z^D%l{k) (32) 

Kpcik, -p, q) D^p{q)D;„{p) T%,{-q,p, -k) , 

where p = k + q. This equation, the Mandelstam equation, is schematically depicted 
in fig. ^. It was already pointed out by Mandelstam that in order to solve this equation 
selfconsistently it is necessary to implement an additional constraint: Without ghosts 
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Figure 3: Diagrammatic representation of the gluon Dyson-Schwinger equation in 
Mandelstam's approximation. 



the Slavnov-Taylor identity (p^) with G = 1 and its solution (25/30) for — > entail 
that 

.t;?„zM=». - ,ymO„-.*(^') = . (33) 

While imposing this additional condition seems consistent with the other assumptions 
in the Mandelstam approximation, it has to be emphasized that concluding (^) as 
a result of (|2^) relies solely on neglecting all ghost contributions in covariant gauges 
(see also the discussion at the end of the next section) . 

The Mandelstam equation in its original form is obtained from eq. (^) upon 
contraction with the transversal projector V^u{k) = 5^^ — k^k^/k"^ and integration of 
the angular variables, 

Z{k^) ^ 1671"^ Jo k^ [Sk^^ 4 A;2 2 J ^'^ ' 
r^'dq^ (7k^ 25 k^ 9\ 2, 

The quadratic ultraviolet divergence in this equation is absorbed by a suitably added 
counter term introduced in order to account for the masslessness condition (|3^), 
see refs. ^, 0. The solution to this equation proposed by Mandelstam's infrared 
analyses proceeds briefly as follows: Assume Z{k'^) ~ l/k"^. On the r.h.s of eq. ( |3^ ) 
this exclusively yields contributions which violate the masslessness condition (|33|). 
Such terms have to be subtracted. Since the kernel on the r.h.s of eq. (^) is linear 
in Z, this is achieved by simply subtracting the corresponding contribution from Z 
in the integrand. Thus defining 

Z{k^) = ^ + C{k'^) , b = const. (35) 



and retaining only the infrared subleading second term in the integrals, eq. ( p4| ) leads 
to a non-linear integral equation for the function C{k^) which, as a solution to this 
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equation, can be shown to vanish in the infrared by some non-integer exponent of 
the momentum [O], 



C(A;2) ~ (A:2)7o ^ = _ i ^ 1.273 , for A:^ ^ . (36) 

Subsequently, an existence proof, a discussion of the singularity structure and an 
asymptotic expansion of the solution generalizing Mandelstam's discussion of the 
leading behavior of C{k'^) in the infrared was given by Atkinson et al. [|^]. 

It was later observed by Brown and Pennington that it is superior for several 
quite general reasons, to be discussed in more detail in the next section, to contract 
eq. ( p2[ ) with the tensor Tl^,^{k) = 5^u — ^k^ky/k"^ . This led to a somewhat modified 
equation for the gluon renormalization function, in particular, without quadratically 
ultraviolet divergent terms P|, 

1 ^ ^ 9' f'^dq^ (7q^ 17 g2 9\ , 



Z(A;2) ^ 167r2 



g2 r^'dq^ hk^ .^M 2^ 
+ tLi i\^-.-^^]Z{<f)- (37) 



167r2 7^2 A;2 

In a previous publication ^ we demonstrated that the solution to this equation has an 
infrared behavior quite similar to the solution of Mandelstam's original equation (pi). 
In particular, explicitly separating the leading infrared contribution according to (|35| ) 
we obtained a unique solution of the form 

Cik^) ~ (A:2)To ^ = ^^229 cos Q arccos " ^ ~ 1-271 , (38) 



for k"^ — > 0. We solved both equations, eq. (|^) as well as eq. (^), using a combination 
of numerical and analytic methods. In the infrared, we applied the asymptotic expan- 
sion technique of ref. Q and calculated successive terms recursively. The asymptotic 
expansions obtained this way were then matched to the iterative numerical solution. 
The results proved independent of the matching point for a sufficiently wide range of 
values, see for details. 

Logarithmic ultraviolet divergences are absorbed in the gluon renormalization 
constant Z3 which can be shown to obey the identity ZgZ^^ = 1 in Mandelstam ap- 
proximation 1^ . This entails that the product of the coupling and the gluon propaga- 
tor, gD^i,{k), does not acquire multiplicative renormalization in this approximation 
scheme. Using a non-perturbative momentum subtraction scheme corresponding to 
the renormalization condition 

Z{e = = 1 (39) 

for some arbitrary renormalization point /i^ > 0, the resulting equation can in both 
cases be cast in a renormalization group invariant form determining the renormal- 
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Figure 4: The gluon renormalisation function gZik'^) = Sna/k'^ + gC{k'^) for eq. (|3 
The dashed Hnes show the corresponding function gC{k'^) for eqs. (|37| ) and (| 



ization group invariant product gZ{k'^) which is equivalent to the running couphng 
g(t,g) of the scheme, 



gZ{k^) = g{tk,g) , 



2 



(40) 



The scahng behavior of the solution near the ultraviolet fixed point is determined by 
the coefficients /3o = 25/2 and 7^ = 25/4 for eq. (|§ vs. /3o = 14 and 7^ = 7 for (|3^) 
which are reasonably close to the perturbative values for Nf = 0, i.e., /?o = 11 and 
7^ = 13/2, the difference being attributed to neglected ghost contributions. 

The full non-perturbative solutions to both equations resulting from Mandel- 
stam's approximation scheme, eq. (34) and eq. ([37|), are compared in figure 0. The 



most important feature of these results to be noted in the present context is their 
implication towards an infrared enhanced quark interaction from the (when ghosts 
are neglected) renormalization group invariant product 

gC{k^) \ 



gD^^ik) = V^y{k) 



Sttct 



+ 



fc2 



(41) 



which would allow to identify the string tension a and relate it to the scale Amom of 
the subtraction scheme Q. The initial question that led to our present study was 
whether the presence of ghosts in Landau gauge could affect these conclusions. 
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It is worth mentioning that this result of the Mandelstam approximation, the 
infrared enhanced gluon propagator (pl|), was in some contradiction with imphca- 
tions from other studies ranging from lattice calculations [ p^ , |35| to implications 
from complete gauge fixings |36|. In addition, an alternative approach to studying 
the Dyson-Schwinger equations of QCD (in Landau gauge) exists which also leads 
to dissentive conlusions |37, 38 1. This approach is based on the observation that an 
exact vertex function and its perturbative estimate differ by expressions which, in 
general, contain essential singularities in the coupling of the same structure as in the 
spontaneous mass scale, see eq. The idea therefore is to improve the perturba- 
tive expansion in approximating the Euclidean proper vertices F by a double series 
in the coupling as well as in the degree of a rational approximation with respect to 
the spontaneous mass scale Aqcd being non-analytic in g'^, i.e., the proper vertices 
are approximated by a truncated double series in functions F^^''^ where r denotes 
the degree of the rational approximation while / represents the loop order of per- 
turbative corrections in calculated from F^'"''^) instead of F^g^,,^ , see ref. |38|. The 
requirement that these non~perturbative terms reproduce themselves in the Dyson- 
Schwinger equations is then used to further restrict the rational Ansatze leaving only 
a few possible solutions. This method preserves the perturbative renormalizibility. 
However, due to the complexity of the resulting equations the approach has to be 
restricted to only a few proper vertex functions. These are then taken as a subset of 
the seven primitively divergent proper vertex functions of QCD: the (inverse) gluon, 
ghost and quark propagators and the 3-gluon, 4-gluon, gluon-ghost and gluon-quark 
vertex functions. Preliminary studies within this approach, restricted to r = 1 and 
Z = for pure QCD in Landau gauge and assuming a perturbative ghost propagator, 
find an infrared vanishing gluon propagator ~ k"^ for small spacelike momenta [^ ]. 
A possible problem of this approach might be the existence of light-cone singulari- 
ties in the 3-gluon vertex function which may be attributed to the low order of the 
employed approximation. A non-trivial ghost selfenergy becomes unavoidable also 
in this approach when perturbative corrections are taken into account via Operator 
Product Expansion techniques [39|. In the meantime several solutions to the order 
r = 1 and / = 1 have been found. However, even though at the present level the 
results might not be regarded as absolutely conclusive, neither of these solutions gives 
rise to an infrared enhancement of the gluon propagator [40U . 

One might think that the particular problem with ghosts in the Landau gauge can 
be avoided using gauges such as the axial gauge, in which there are no ghosts in the 
first place. As for studies of the gluon Dyson-Schwinger equation in the axial gauge 



in, 11, O] it is important to note that, so far, these rely on an assumption on the 



tensor structure of the gluon propagator. In particular, an independent additional 
term in the structure of the gluon propagator has not been included in presently 
available studies of the gluon DSE in axial gauge |1S]. This term vanishes in per- 
turbation theory. If, however, the complete tensor structure of the gluon propagator 
in axial gauge is taken into account properly, one arrives at a coupled system of 
equations which is of considerably higher complexity than the ghost-gluon system in 
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the Landau gauge |14|. Some preliminary progress in this direction has so far only 
been obtained in light-cone gauge |41]. Available studies of the gluon propagator 
in the axial gauge can therefore not be regarded any more conclusive than those in 
Landau gauge. In fact, the present situation in axial gauge DSE studies can be sum- 
marized to be quite comparable to the Mandelstam approximation in Landau gauge: 
Early truncation and approximation schemes to the gluon DSE in axial gauge, due 
to the simplified tensor structur, lead to equations quite similar to Mandelstam's 
equation |0, |ll|. Furthermore, these studies led to an infrared enhanced gluon 
propagator, D{k) oc a/k^ for k'^ — > 0, analogous to the solution to Mandelstam's 
equation |^ Subsequent studies came to a somewhat dissentive conclusion 

obtaining an infrared vanishing gluon propagator in the identical scheme in axial 
gaug e p!^ . This was, however, later attributed to a sign error in the resulting equa- 
tion |13|. Therefore, previous DSE studies of the gluon propagator in Landau gauge 
as well as in the conventional axial gauge seemed to agree in indicating an infrared 
enhanced gluon propagator oc fj/Zc^. On one hand, such a solution is known to lead 
to an area law in the Wilson loop I^H, on the other hand it was also argued from 
positivity constraints that the full axial gauge gluon propagator should not be more 
singular than l/Zc^ in the infrared [^]. While this apparent puzzle might be resolved 
by using suitable principle value prescriptions for the axial gauge singularities 1 44 ] , it 
is to be expected that the complete tensor structure of the gluon propagator in the 
axial gauge will have to be fully taken into account before a more conclusive picture 
can emerge. As for the Landau gauge, we will find that, instead of the gluon prop- 
agator, the previously neglected ghost propagator assumes an infrared enhancement 
similar to what was then obtained for the gluon. 

Progress is desirable in axial gauge as well, of course. A further prerequisite 
necessary for this, however, will be a proper treatment of the spurious infrared diver- 
gences which are well known to be present in axial gauge due to the zero modes of 
the covariant derivative |44]. This can be achieved by either introducing redundant 
degrees of freedom, i.e. ghosts, also in this gauge |Q (by which it obviously looses its 
particular advantage) or by using a modified axial gauge |^^, which is specially de- 
signed to account for those zero modes. Ultimately, progress in more than one gauge 
will be the only reliable way to asses the influence of spurious gauge dependencies. 



4 Infrared Dominance of Ghost Contributions 

From now on we concentrate on the coupled system of integral equations derived in 
section 2 for the propagators of gluons and ghosts. Instead of attempting its direct 
numerical solution, we will solve this system simultaneously in an one-dimensional 
approximation which allows for a thorough analytic discussion of the asymptotic in- 
frared behavior of its solutions. This is necessary in order to obtain stable numerical 
results from a matching procedure which generalizes the techniques successfully ap- 
plied to the gluon DSE in Mandelstam approximation The approximation we use 
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to arrive at the one-dimensional system of equations is thereby designed to preserve 
the leading order of the integrands in the infrared limit of integration momenta. At 
the same time, it will account for the correct short distance behavior of the solutions 
(behavior at high integration momenta). Because of the particular simplicity of the 
ghost DSE, we start with this equation to introduce and motivate a modified angle 
approximation. From (^) with the vertex ( p2|) we obtain the following equation for 
the ghost renormalization function G{k'^), 

' ~ ^^nJ ^Am,)^-^ (42) 



Here and in the following we use that in Landau gauge Zi = 1 |28|. In order to be 
able to perform the integration over the 4-dimensional angular variables analytically, 
we make the following approximations: 

For q^ < A;2 we use the angle approximation for the arguments of the functions 
Z and G, i.e., G{p^) = G{{k - qf) G{k'^) and Z{p'^) ^(A:^). This obviously 
preserves the limit q^^ ^ q of the integrand. 

For g2 > /j2 make a slightly different assumption, which is that the functions 
Z and G are slowly varying with their arguments, and we are thus allowed to replace 
GijP') ~ G{k'^) G{q'^). This assumption ensures the correct leading ultraviolet 
behavior of the equation according to the resummed perturbative result at one-loop 
level. For all momenta being large, i.e. in the perturbative limit, this approximation 
is well justified by the slow logarithmic momentum dependence of the perturbative 
renormalization functions for ghosts and gluons. Our solutions will resemble this 
behavior, justifying the validity of the approximation in this limit. 

With this approximation, we obtain from (^) upon angular integration. 



Gik"^) 167r2 4 [Jo k"^ k"^ ^ ' ^ ' Jk^ q^ 

where we introduced an 0(4)-invariant momentum cutoff A to account for the loga- 
rithmic ultraviolet divergence, which will have to be absorbed by the renormalization 
constant. 

A preliminary discussion of the implications of eq. ( |4^ ) on the infrared behavior 
of the renormalization functions reveals the following: 

Making the Ansatz that for re := fc^ — > the product of G and Z behaves as 
Z{x)G{x) ~ x'^ for K 7^ 0, it follows readily that 

G{x) ~ x'"" and Z{x) ~ x^'' as x^O . (44) 
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Furthermore, in order to obtain a positive definite function G{x) for positive x from 
a positive definite Z{x), as x — > 0, we find the necessary condition 1/k — 1/2 > 
which is equivalent to 

< K < 2 . (45) 



The special case k = leads to a logarithmic singularity in eq. ( [43D for x ^ 0. 
In particular, assuming that ZG = c with some constant c > and x < xq for a 
sufficiently small xq, we obtain G^^{x) — > c{3Ncg'^ /64tt'^) ln{x/xQ) + const and thus 
G{x) — > 0~ for X ^ 0, showing that no positive definite solution can be found in this 
case either. 

The gluon DSE (^) is more complicated. In a manifestly gauge invariant formula- 
tion the gluon Dyson-Schwinger equation in the covariant gauge would be transverse 
without further adjustments. For the following decomposition of all its contributions 
to the inverse gluon propagator: 

D-^{k) = A{k^) 6^, - B{k^) k^k, + 1^ k^k, , (46) 

this implies that A{k'^) = B{k'^) = Z~^{k'^). The longitudinal part of the gluon prop- 
agator does not acquire dressing and cancels with the one of the tree-level propagator 
on the r.h.s. of the gluon DSE (^) (^ = Z^^^q)- Although we can do our best to use 
implications of gauge invariance in constructing vertex functions, we cannot expect 
to arrive at an exactly gauge covariant truncation scheme. This fact is reflected in 
^ 7^ -B as obtained from truncated gluon DSEs. An additional source of spurious 
longitudinal terms in the gluon DSE is the regularization by an 0(4)-invariant Eu- 
clidean cutoff A which violates the residual local invariance, i.e., the invariance under 
transformations generated by harmonic gauge functions {d'^A{x) = 0). The straight- 
forward elimination of spurious longitudinal terms by contracting eq. (|2|) with the 
transversal projector Vfj,y{k) is known to result in quadratically ultraviolet divergent 
contributions which are of course artifacts of the regularization not being gauge in- 
variant. As observed by Brown and Pennington [^], in general, quadratic ultraviolet 
divergences can occur only in A{k'^). Therefore, this part cannot be unambiguously 
determined, it depends on the momentum routing. An unambiguous procedure is to 
isolate B{k'^) by contracting the truncated gluon DSE with the projector 

n^u{k) = S^,-A'^ , (47) 

and to set Z{k'^) = 1/B{k'^). Note that an additional advantage of this prescription 
is that the constant contribution of the tadpole term in the full gluon DSE does not 
enter in B{k'^) either. Since the non-perturbative tadpole does not necessarily vanish, 
it is in fact an example of a contribution which, if neglected, can lead io B. From 
these various reasons it is clear that one should concentrate on B rather than A in 
truncated gluon DSEs [||, g. 



20 



We gave the detailed justification for this, because it imphes that the usual ar- 
gument of the irrelevance of longitudinal terms in the ghost-gluon vertex in Landau 
gauge does not apply to the ghost-loop in the gluon DSE (|2|). We are interested in 
just the unambiguous term B proportional to fe^fcjy. Therefore, contracting (§) with 
TZ{k), we see that we have to use the full form of the vertex given in eq. (^) in the 
gluon DSE. With this and the 3~gluon vertex as given in eqs. (|25|/ l26|) we obtain. 



- Zi—— / ——- <Ni{p ,q ■,k ' 



Z{k^) ^ 6 y (27r)4 [ ' ^ Z{k^)G^{k^) 

...2 2 ,2^Z{p^)G{p^) ,,,2 2 ,2^Z{q^)G{q'^)\ G{k^ 



+ 



G(g2) ' / (^(^2) j ^2_p2^2 



(2vr) 



4 



[qn{k)q)[G{k')G{p') - G{q^)G{p')) 
-[qn{k)p)G{k')G{q')\^^ . (48) 



The functions Ni{x,y;z) = Ni{y,x;z) and N2{x,y;z) are given in appendix for 
completeness. These functions have singularities for coinciding arguments canceling 
only in their sum, Ni{x, y; z)+N2{x, y; z) + N2{y, x; z), which determines the 3-gluon 
loop contribution in the Mandelstam approximation |§, |^ |8|, |9| . In the present case, 
the modified angle approximation allows us to combine the three terms of the 3-gluon 
loop for > k'^ in the same way as in Mandelstam approximation. For q^ < k'^, 
however, the second term in the 3-gluon loop in (|4^) can only be combined with the 
other two by encountering an additional error of the form 

Zo'^'^ f ^''^ 7Vr„2 ,2.,2^/^(pW) Ziq^)G{q') \ G{k^) 

In this contribution the singularity in N2{x, y; z) ~ l/{x — y) is cancelled by the terms 
in brackets which vanish for q'^ p^. However, employing an angle approximation on 
these terms, the cancellation of the singularity is destroyed and an artificial singular 
contribution would arise. This demonstrates that approximations have to be used 
with some care in order to avoid spurious divergences (and imaginary parts). The 
angle approximation cannot be reasonably applied to the contribution (^9|). We will 
omit this contribution to the gluon DSE for q^ < k"^ in addition to the modified angle 
approximation, since its inclusion would preclude a one-dimensional reduction of the 
equations. 

Furthermore, this additional approximation has no influence on the motivation 
for the angle approximation for momenta q^ < k"^ , i.e., to preserve the infrared limit 
of the integrands, because all contributions of the 3~gluon loop to the r.h.s in eq. 
(^) are subleading in the infrared as compared to the contributions of the ghost 
loop, as our preliminary discussion of the infrared behavior of the solutions reveals 
(see below). We will therefore put particular emphasis on the ghost contributions to 
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the coupled system of Dyson-Schwinger equations. In order to isolate their effect, 
it seems reasonable to treat the 3-gluon loop as analogous as possible to previous 
studies of the gluon DSE in the Mandelstam approximation. This is achieved by the 
modified angle approximation to (^8|), if in addition, for < k"^, the contribution 
(E^) in the 3-gluon loop is omitted. This yields from (E^ upon angular integration. 



^3 + ^iT7T-o^W -i:T{ni:i-^i:9.-^]Z(<i)G{q 




167r2 3 lio A;2 2 A;2 



,ri^(if-^)^toW)) (50) 

3 2 Jk'i 



The only difference in the 3-gluon loop as obtained here versus the Mandelstam 
approximation is that the gluon renormalization function Z is replaced by the prod- 
uct ZG, see eq. ( p7[ ) in sec. 3. This is one consequence the presence of the ghost 
renormalization function has on the dressing of the 3-gluon vertex. The system of 
equations ([43| ) and (^) is a direct extension to the gluon DSE in the Mandelstam 
approximation. As such it is very well suited to investigate the influence of ghosts 
on conclusions of previous studies based on this approximation. 

To asses the relative importance of the neglected term (49) we calculated this con- 



tribution without angle approximation using, however, the selfconsistent solutions for 
Z(A;2) and G{k'^) as obtained from the one-dimensional set of equations ( ^3[ ) and ([50|), 
c./., the numerical method and results given in section 8 below. In figure ^ we plot 
the resulting contribution in comparison with the terms retained on the r.h.s. of eq. 
(^) which, for selfconsistent solutions, are equivalent to 1/Z{k'^). Though even small 
terms can, in principle, have a considerable effect on the non-linear selfconsistency 
problem, the fact that the additional contribution to the gluon DSE due to ( ^9D is 
comparably small for all momenta supports this additional approximation. In partic- 
ular, neither the infrared nor the ultraviolet asymptotic regimes are affected by this 
contribution at all. 

Before we turn to the renormalization of equations ( ^3| ) and ( |50| ) in the next 
sections, we conclude this section with an extension of our preliminary infrared dis- 
cussion. From (E^ we find that ZG cx'^ for x ^ implies. 



^ (51) 

Z{x) ^ (f'^o" - f )) cV" , 0<K<2, (52) 

where is the leading order perturbative coefficient of the anomalous dimension 
of the ghost field. Accordingly, the ghost-loop gives infrared singular contributions 
~ x~2'^ to the gluon equation (pO) while the 3-gluon loop yields terms proportional 
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Figure 5: The dismissed contribution ( ^9|) (dashed hne) compared to the contribu- 
tions retained on the r.h.s of (|50|), i.e., XjZiy?) (sohd Hne) as a function of x = I? ja 
with a ~ (350MeV)2, c./., sec. 8. 



to x*^ as X — > 0, which are thus subleading contributions to the gluon equation in the 
infrared. With eq. (|5l| ) the leading asymptotic behavior of eq. (5C) for x — > shows 
that 
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0.92 . 



(53) 



(54) 



In the last line we compared (53) to ( |52[ ) and used the restriction (|45) obtained from 
the ghost equation on the exponent k (ruling out the other sign in (5J)). Note that 
K is independent of the number of colors A'^c- 

This leading behavior of the gluon and ghost renormalization functions and thus 
their propagators is entirely due to ghost contributions. The details of the treatment 
of the 3-gluon loop have no influence on above considerations. This is in remarkable 
contrast to the Mandelstam approximation, in which the 3-gluon loop alone deter- 
mines the infrared behavior of the gluon propagator and the running coupling in 
Landau gauge 0, ^, ^J. As a result of this, the running coupling as obtained from 
the Mandelstam approximation is singular in the infrared In contrast, as we will 
show in the next section, the infrared behavior derived from the present truncation 
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scheme implies an infrared stable fixed point. This is certainly a counter example 
to the frequently quoted assertion that the presence of ghosts in Landau gauge may 
have negligible influence on physical observables at hadronic energy scales. 

With the infrared behavior of gluon and ghost propagators obtained above, a 
comment on its implications for the vertex functions is in order. Starting with the 



ghost-gluon vertex (21) we realize that the limit of vanishing ghost momenta is 
regular: 

G^{p,q)^ -ik^ for p^O (55) 

Gf,{p,q)^ for g^O , (56) 

where we used that G{k^) ~ {k'^)^'^ for fc^ — > 0. On the other hand, for vanishing 
gluon momentum A: — > the vertex diverges as 

G,(,,<,)-2ip,§||~JPl^ lor .-0 . (57) 



Since the exponent k is a (positive) irrational number, see eq. (p4[), the corresponding 
divergence cannot be interpreted to reflect the presence of a massless excitation. This 
divergence is in fact weaker than a massless particle pole (k < 1) and presumably 
lacking a physical interpretation. 



Similarly, the 3-gluon vertex as given in equations (25/^) shows analogous in- 
frared divergences resulting from the infrared enhanced ghost renormalization func- 
tion. 



^tj.up{p, q, k) G{k'^) < {ipp5up + ipu^tip - 2ipp5^i^ }^ (58) 




+ 2^Ppp2p,.(p)(^pj|^ + ^5^)^ for fc^O 

This is the correct result to satisfy the differential Slavnov- Taylor identity obtained 
from eq. ( ]2^ ) for /c — > 0. Note that in the previous studies of the gluon DSE in 
Mandelstam approximation |6|, 1^, |8|, , the requirement for the 3-gluon vertex func- 
tion to obey the differential Slavnov-Taylor identity led to the so-called masslessness 
condition, eq. (|33|), p^ /Z{p'^) — > for p^ — > 0. In absence of ghost contributions, the 
gluon DSE had to be supplemented by this as an additional constraint. This original 
condition is violated by the infrared behavior of the gluon propagator found here, 
Z{x) x^'^ . The correct replacement of this condition for the present case using 
(P^) with the solution (^/^6|) is, however. 



pW)-0 and ^^-^^^^0 for _ . (59) 

These two necessary conditions are obeyed by the infrared behavior obtained here, 
i.e., G{x) ~ 1/ {G{x)Z{x)) ~ x~'^ with n ~ 0.92 for a; — > 0, without further ad- 
justments. We see, however, that generally the original bounds on k obtained from 
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the consistency of the ghost DSE, when combined with conditions (^), have to be 
restricted further to < k < 1. The possibihty that the masslessness condition in its 
original form might not hold for non-Abelian gauge theories has been pointed out 



previously 48]. Our present results support this conjecture. 

Finally, since we do not expect quarks to screen the divergence in the ghost renor- 
malization function completely, its presence in the quark-gluon vertex function, see 
eq. (p9|), will have a similar effect there, too. We will demonstrate on the example of 
the running coupling in the next section, how these apparently unphysical divergences 
in the elementary correlation functions of gluons, ghosts and quarks can, in princi- 
ple, nevertheless cancel in physical quantities. The lack of a physical interpretation of 
these divergences (other than maybe reflecting confinement) should, however, not be 
too surprising for a Euclidean field theory violating reflection positivity (see sec. 11). 
Note also that in all skeletons of the kernels in relativistic bound state equations for 
hadrons, the combination of dressed quark-gluon vertices with the non-perturbative 



gluon propagator will give contributions of the form (c./., eqs. (|28|j2 

~ G\k^) D^,{k) Tl%p,p -k)® Tl%q -k,q), (60) 

where r^^(p, q) stands for a vertex function with non-perturbative dressing of similar 
structure to what occurs in an Abelian theory, e.g., a Ball-Chiu construction from its 
Ward-Takahashi identity. This combination of the ghost renormalization functions 
with the gluon propagator, g'^G'^{k^)D{k), in the effective quark interaction is free 
of the unphysical infrared divergences for k ^ 0. It is furthermore independent of 
the renormalization scale in Landau gauge and suited to define a non-perturbative 
running coupling as we discuss next. 



5 Subtraction Scheme and Non— Perturbative Running 
Coupling 

Before we discuss the renormalization of the Dyson-Schwinger equations for glu- 
ons and ghosts in the next section, some introductory remarks on the choice of the 
non-perturbative subtraction scheme and its relation to the definition of the run- 
ning coupling are necessary. In particular, we will see that the preliminary infrared 
discussion of the last section already yields an important first result: it implies the 
existence of an infrared fixed point. 

We begin with the identity for the renormalization constants (as introduced in 
eq. (D), 

Zi = ZgZl^^Zs = 1 , (61) 

which has been shown to hold in the Landau gauge [^8|. In fact, our Slavnov-Taylor 
identity (|T^ ) proves that it remains valid in general covariant gauges to the extend 
that irreducible 4-ghost correlations are neglected. Prom eq. (^Tj) it follows that 
the product g'^Z{k^)G'^{k'^) is renormalization group invariant. In absence of any 
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dimensionful parameter this (dimensionless) product is therefore a function of the 
running couphng g, 

g^Z{k^)G\k') = f{f{tk,g)), tj, = Une/i?. (62) 

Here, the running couphng g{t,g) is the solution of d/dtg{t,g) = I3{g) with g{0,g) = 
g and the Callan-Symanzik /3-function (3{g) = —(3og^ + 0{g^). The perturbative 
momentum subtraction scheme is asymptotically defined by /(x) x for j; — > 0. 
This is realized by independently setting 

Z{n^) = 1 and G(^2) = 1 (63) 

for some asymptotically large subtraction point k'^ = ji'^. If the invariant product 
g"^ Z {k'^)G'^ {k'^) is to have a physical meaning, e.g., in terms of a potential between 
static color sources, it should be independent under changes (g,/x) — > (g',/i') accord- 
ing to the renormalization group for arbitrary scales n'. Therefore, 

g'Z{f^'')G\fi'') = g'^ = 52(ln(^7^),5) , (64) 

and, fix) = X, Vx. We adopt this as a physically sensible definition of a non- 
perturbative running coupling in the Landau gauge. 

This definition is an extension to the one we used in the Mandelstam approxima- 
tion, eq. (^0|). In ref. we proved in this approximation without ghosts the identity 
ZgZ-^ = 1 implying that gZ{k'^) is the renormalization group invariant product in 
this case. The according identification of this product with the running coupling 
gZ{k'^) = g{tk,g) is equivalent to the non-perturbative renormalization condition 
Z{ij?) = 1,V;U in Mandelstam approximation (in which there is no ghost renormal- 
ization function). 

In the present case, it is not possible to realize f{x) = x by simply extending the 
perturbative subtraction scheme (]63| ) to arbitrary values of the scale ^, as this would 
imply a relation between the functions Z and G which is inconsistent with the leading 
infrared behavior of the solutions discussed in the last section. For two independent 
functions the condition ( |63|) is in general too restrictive to be used for arbitrary 
subtraction points. Rather, in extending the perturbative subtraction scheme, one is 
allowed to introduce functions of the coupling such that 

Z{li^) = fA{g) and Gifi^) = fcig) with /^/^ = 1 , (65) 

and the limits Ja.g ~^ 1 > 0. Using this it is straightforward to see that for 

A;2 / ^2 Qj^g j^j^g ^ (lnfeVA'^)/2), 

Z{k^) = exp I - 2 1'^'"'^ dl ^1 fAmk.g)) , (66) 

Gik") = exp|-2^'^*"'^d/^|/G(5(tfc,<7)). 



26 



Here 7^1 ((7) and 'ycid) are the anomalous dimensions of gluons and ghosts respectively, 
and P{g) is the Callan-Symanzik /3-function. Eq. ( |6T|) corresponds to the following 
identity for these scaling functions in Landau gauge: 

2ig{9) + 1a{9) = -- m ■ (67) 
9 

Therefore, we verify that the product g^ZG^ indeed gives the running coupling {i.e., 
eq. ( |62D with f{x) = x). Perturbatively, at one-loop level eq. ( p^ ) is realized 
separately, i.e., jg{9) = -^P{9)/9 and 7^(5) = -(1 - '^^)P{9)/9 with 5 = 9/44 for 
Nf = and arbitrary Nc. Non-perturbatively we can separate these contributions 
from the anomalous dimensions by introducing an unknown function e{g), 

jg{9) =:-iS + e{g)) ^ 7^(5) = "(1 -25- 2eig)) ^ . (68) 



This allows us to rewrite (|6^) as follows: 

= (l!(^)-^.p{_4/;""^,f}/,(5fe..)), (69) 

This is generally possible, also in the presence of quarks, in which case 6 = 7o'//5o = 
9A'^E/(44A'^r; — 8Nf) for Nj flavors in Landau gauge. The above representation of the 
renormalization functions expresses clearly that regardless of possible contributions 
from the unknown function €{g), the resulting exponentials cancel in the product 
G^Z. For a parameterization of the renormalization functions, these exponentials 
can of course be absorbed by a redefinition of the functions Ja.g- The only effect 
of such a redefinition is that the originally scale independent functions fA,G{9{tk,9)) 
will acquire a scale dependence by this, if e / 0. We therefore use the form ( |69[ ) to 
motivate the following parameterization for G and Z, 

Z{k') = (70) 
f F{x)\^ 1 

G{k ) = I j with X := k'^/a and s := fi'^/cr , 

where a is some currently unfixed (RG invariant) scale parameter. This obviously 
implements the renormalization condition G^Z'|^2=^2 = 1 and from the definition 
of the running coupling (^) we find that (p'{tk,g) ~ We fix the constant of 

proportionality for later convenience by setting (with I3q = llA''c/(487r^) for Nj = 
quark flavor), 

q2 1 

Pof{tk,9) = F{x) and as{^i) = ^ = ^(*) ' (^1) 
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accordingly. As stated above, in general, a non-trivial contribution from e{g) will re- 
sult in the function R being scale dependent. For the present truncation scheme it will 
be possible, however, to obtain explicitly scale independent equations to determine 
the running coupling ~ F as well as the function R thus showing that the solutions 
for the renormalization functions G and Z obey one-loop scaling at all scales. In 
particular, it implies that the products ^^^G and g'^^^~'^^^ Z are separately renormal- 
ization group invariants in the present scheme (as they are at one-loop level). As for 
the renormalization scale dependence, the non-perturbative nature of the result is 
therefore buried entirely in the result for the running coupling. 

We are now in a position to discuss the implications of the preliminary results for 
the infrared behavior of the solutions G and Z at the end of the last section without 
actually solving equations (50,E^). From eqs. (|5l|) and ([5^) we find for k"^ — > 0, 



g'Zik^)G\e) = -9\t,,g) '^^^ [j^ (i - i))"' =: gl . (72) 



The critical coupling scales with the number of colors as g^ ~ 1/-/Vc implying that 



in our approach gf.Nc is constant. This agrees with the general considerations of the 
large A'^c-limit in which g'^Nc is kept fixed as Nc becomes large |49]. For Nc = ?> and 



with eq. ( |54D for k we obtain g^ ~ 119.1 which corresponds to a critical coupling 
Oc = 5^/(47r) ~ 9.48. This is a remarkable result in its own, if compared to the 
running coupling as it was analogously obtained from the Mandelstam approximation 
Q . The dynamical inclusion of ghosts changes the infrared singular coupling of the 
Mandelstam approximation to an infrared finite one implying the existence of an 
infrared stable fixed point. 



6 Renormalization 

We discuss the actual renormalization of equations (^) and (^) once again beginning 
with the simpler and thus more evident case of the ghost DSE. To see the effect of 
the renormalization condition we use the parameterization ( [70| ) of the last section in 
the ghost DSE (||) and set k"^ = fi"^ x = s) with Pog"^ = F{s) and obtain, 

R{s) = Zs-6 (^lRis)Fis) + F\s) ^ R{y)F^-\y)^ , (73) 
where L := h? /a \s the ultraviolet cutoff and S = / l3o ( = 9/44 for iVj = 0). This 



determines the renormalization constant Z3 for arbitrary scales s. We can use (73) 
to eliminate Z^ in (Bsh and obtain analogously for general x ^ s, 



Rix) Ris) ^^^^ 



F^{x) F^{s) 

5i^^{R{s)F'-'{s)-R{x)F'-'{x)) + J^' ^ R{y)F'-\y)^ 
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This equation is now ultraviolet as well as infrared finite. Note that the preliminary 
infrared analyses corresponds to 

F{x) -^a:= Po9c- 8.3 and R{x) -^bx"" , x^O , (75) 

where b is some currently unfixed constant depending on the choice of the scale 
parameter a which we did not specify so far. Equation ( [7^ ) is valid for all scales s, 
in particular, we may let s ^ to obtain, 

^ = ' (r 7 ^(y)^^-^(^) - I Ri-)F''H^)) ■ (76) 
We can use this equation in ([f3|) exchanging x <-> s to write. 



Zs = F'{s) (^6 ^ i?(y)Fi-'^(y) j . (77) 

If we use this form for in (|4^), we arrive directly at (^) which shows that it 
remains valid also for non-zero scales s (we have thus shown that all x-dependent 



terms cancel separately from the s-dependent ones in eq. (74) with the integration 
range (s, x) split into (0, s) and (0, x)). It is instructive to replace the integral in ( |77| ) 
by again using eq. ( |76D now with x L, 

«<^' f,^inL;)F%,)^^-f:^JS^]\ (78) 



for the perturbative limits of large x, F{x) — > 1/lnx and R{x) — > 1, showing nicely 
the scaling limit, in particular, the possible interpretation of the bare coupling as the 
running coupling at the cutoff scale in the present scheme. Furthermore, the renor- 
malization constant in the scaling limit (large L) becomes identical to the multiplica- 
tive factor for finite renormalization group transformations of the ghost propagator, 
Z^{^'^ , pi'"^) for /i — > /i', at the cutoff scale, Z3 = Z^{^'^,L). This detailed study of the 
scaling behavior of the ghost DSE justifies the use of the improved angle approxima- 
tion for integration momenta > k"^ in this case. We now turn to the gluon DSE 
using the discussion above as a guideline for the anticipated scaling behavior of the 



gluon equation (pOl). 

The gluon DSE has the additional problem of the presence of the renormalization 
constant for the 3~gluon vertex Zi in the gluon loop. We will address this issue in 
the following. The first step is to rewrite (50) with ( |70[) and set x = s, 

3 F^s) 
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where we have used (77) to express the logarithmically divergent term in the 3- 
gluon loop in terms of the constant Z3. Furthermore, we used that the integral 
Sx{dy/y'^)RF^-^ J^{dy/y^)RF''-^ is ultraviolet finite for functions R and F 
resembling the perturbative behavior at high momenta. We now have to specify 
Zi. We will discuss three possibilities in the following. The first and probably most 
obvious will lead to a contradiction. The numerical results of the other two will 
not show major qualitative differences except at short distances. Comparing to the 
known perturbative behavior at short distances, we suggest that the third method of 
treating Zi is suited to restore some of the effects of the truncation which are at the 
root of the present problem. 

a) Zi = Zs/Zs 

The most rigorous way would be to use the Slavnov-Taylor identity Zi = Z^/Z^ 
which follows from Zi = 1 in Landau gauge. We could then attempt to eliminate Z3 
from (|5^) in analogy to ( [74[ ) for x / s writing, 

1 1 



ZsF'jL) 1 [ r-dy (7_y^ 
^1-2-5(5) 11 |io X \2x^ 

+ xl r^R{y)F'-\y) 
o Jx y 




(80) 

R{y)F^-'{y) 




F\x) dy y F^y) F\s) dy y F\y) 



R{x) Jo XX R{y) R{s) Jo s s R{y) 
i(F^\x) F2'5(s)\ 1 dy F^\y)' 



3 \ R'^{x) R-^{s) J 2js y R'^{y) j ' 

where we used F^{s)Zi = Z^F^{L) from (^) for large L. This factor ZsF^L) is the 
only cutoff dependence left in (|80|). Expecting a scaling limit as in the ghost case, i.e., 
Z3 — > {F{s)/F{L))^~^^ , we see that the prefactor of the 3-gluon loop contribution 
to (|0|) behaves as ~ 1/F^-^^{L) which is singular in the limit L— >oo, ifl — 35>0. 
Perturbatively we expect this to be the case since 36 = Sjq / Po < 1 <^ /3o > 27/4 
which is true for A'^^ < 7. There is no term left in the otherwise ultraviolet finite eq. 
(pO|) to absorb the singularity in this case. 

It is interesting to note that perturbatively for Nf > 7 the leading anomalous 
dimensions are such that Zi = Z3/Z3 — > in the scaling limit in Landau gauge. 
This implies that the contribution of the 3-gluon loop to the gluon DSE vanishes 
completely for Nf > 7. 

In order to reproduce the leading logarithmic behavior from perturbation theory 
for high momenta, being determined by 1 — 3(5 = 17/44 for Nf = 0, we find that it is 
not possible to renormalize the DSEs for gluons and ghosts in the present truncation 
scheme while retaining the identity Zi = Z^/Z^. Could it alternatively be possible to 
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find a consistent solution for the special case 1 — 35 = 0? Obviously, eq. pQ ) would 
have a finite scaling limit in this case. However, a detailed study of the asymptotic 
behavior of eqs. ( [f6[ ) and ( pO| ) shows that this leads to a contradiction. From eq. {7t) 
for large x with i? ^ 1 we obtain a differential equation for F asymptotically, which 
has a solution F{x) = 1/ {const. + lnx) with 5 = 9/(4/3o). Therefore, in order to have 
5 = 1/3 we need /3o = 27/4. The same sort of solution is possible asymptotically from 
(^), however, this time we find 5 = 13/(2/?o) = 1/3 which gives the contradiction. 



The coefficients of the ultraviolet dominant terms of eqs. (76) and ( pO| ) are consistent 
only with the perturbative values for 70", and (3^. We therefore conclude that it 
is necessary to stay with the perturbative short distance character of the equations 
and abandon the identity Zi = Z^/ instead. These considerations show that our 
truncation scheme is unable to obey the implications of gauge invariance exactly 
which may not be surprising after all. However, we will see in c) below that the 
hope to obey all Slavnov-Taylor identities for the renormalization constants was 
only slightly too optimistic. 



h) Zi = l 

The rigorous case above being inconsistent for the present truncation scheme, it 
may seem natural to try the easiest next, setting Zi = 1 since this constant is not 
necessary to absorb any ultraviolet divergences at the present level anyway. We will 
see below that Zi = 1 \s essentially equivalent to setting it to an arbitrary (finite) 
constant cz^- As in the case of the ghost DSE, we look at the limit of eq. ( [79| ) for 
s — > in which this equation is singular in contrast to the ghost eq. (|7^). Using the 
limits ( [75| ) the leading behavior of eq. ( [79| ) up to contributions vanishing for s ^ 
is given by. 
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From the definitions of k and a = /3o5c ' ^"^^ 



^ dy F^\y) 
y R\y) 

we find 



(81) 



a 
11 

we obtain 
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a 1 1-1 

and with 
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(82) 



We assumed that Z3 = Z^[s) is finite in the limit s ^ which is necessary from eq. 
(|6l]), since due to the infrared fixed point Zg and Z^ (see eq. (|77|)) are infrared finite. 



Furthermore, the above integral exists only if F^° jR^ a /{b y ) + 0{y'^) , r > 0, 
which we will verify when we come to discuss the infrared expansion of the solutions 
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for R and F. We are now able to eliminate all ultraviolet divergences from the gluon 
DSE which can be written as 



I { r dy h y"^ 17 y 9 x\ \pi~5f 
i?2(x)Fi-2'5(x) ~ ^^U\Jo~[2l^~Yx~8^ y]^^^ 



^ f°° dy ^, A 1 \2,F%x) rdyyFHy) 

l F'^\x) _ I dy f F^\y) _ a^^ ^ ^ ^ 



3 R'^{x) 2 Jo y \R'^{y) b^y^"" J 4k b^x'^'' j 

1 /3 

We see from this equation that with a redefinition iZ' := c^^ R and accordingly 
1 /3 

6' = c^-i^ 6 we can eliminate the finite constant cz^ from (p^). The ghost DSE ( [7^ ) is 
linear in R and thus in its form unchanged by this redefinition. It implies a scaling 
of the propagators by a multiplicative constant which does not appear in physical 
quantities like the running coupling and is thus unessential. We will therefore solve 
the coupled system ( |76| ) and (^) for Zi = cz^ = 1. It is easy to see from those terms 
in eq. ( |83|) which are dominant for large x, that the resulting gluon propagator cannot 
be expected to have the correct perturbative ultraviolet behavior (its perturbative 
anomalous dimension): The leading ultraviolet contribution from the ghost loop has 
the wrong sign and should be compensated by the according contribution from the 
3-gluon loop. However, in eq. (|8^) the latter gives a contribution qualitatively equal 
to the asymptotic behavior of the ghost DSE ( [76|) which is logarithmically subleading 
as compared to the ghost loop (since Z3/Z3 for Nf < 7). It is clear that we 
cannot expect the solutions to resemble the perturbative behavior at asymptotically 
high momenta for this reason. 

Evidently, the insufficiencies of the truncation scheme appear here in the approx- 
imations used to simplify the 3-gluon loop. We reiterate, however, that the 3-gluon 
loop does not affect the leading behavior of the ghost-gluon DSEs in the infrared. 
We suggest a further manipulation in the following which is motivated by the desired 
short distance behavior of the solutions and which is suited for its restoration. We 
will compare our numerical results obtained for Z\ = \ \,o those from the procedure 
outlined in the following. In particular, this will verify that the additional manipu- 
lation improves the short distance (high momentum) behavior of the solutions, and 
that this is indeed all it does (not affecting conclusions on the infrared behavior). 

c) Zi = (F(.)/F(y))i-3^ 

As a first remark, we note that we found Z3 = (F{s)lF(L)f in the scaling 
limit. Analogously, we would like to obtain Z3 = {F(s)/F{L)y^'^^ and Zi = 
{F(s)/F{L))^~^^ from the corresponding Slavnov-Taylor identity. We saw above 
that the latter, i.e., Z\ = Z3/Z3, was not consistently possible at the present level 
of truncations. Therefore, we suggest to set Zi = {F{s)/F{y))^~^^ instead, where 
y = q'^/a denotes the loop momentum. We repeat the same steps as before to elimi- 
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nate the renormalization constants and, in place of eq. (pO[), we arrive at 
1 1 



i?2(a;)Fi-2'5(x) R^{s)F^-^^{s) ^^^^ 

11 yjo X yix"^ 2x8 y J o Jx y 

— (x ■f-^' s)| + ghost loop contributions , 

where we did not explicitly repeat the unchanged contributions from the ghost loop 
again (see eq. (pO|)). If we now repeat the exact same steps as in the case of the ghost 
DSE, the limit s ^ now leads to 

11 rdyf7y^ 17 y 9 , or, , 

— - ^- - o + 7- R{y)F^^{y) 



i?2(a;)Fi-25(2;) Jo X \2x^ 2x8 y^ 

3 F\x) dy y F^y) 1 F'^\x) 




y2 2 R{x) Jo X X R{y) 3 R'^{x) 

4k iP'x'^'^ 



(85) 



which can again be use to obtain the gluon renormalization constant, exchanging 
variables x ^ s, with the result, 

7 pi-^r ^l^ f'^'^y T,( , I ["-dy ( F^\y) a'' \\ 

^'^"^ ^'\nJo ^""^'^^ ^^)-^io 7iw-6V^jj- ^''^ 



This is analogous to eq. ( |78| ) and valid for all s. Using eq. (|8^) with x ^ L for large 
-L we furthermore obtain the desired scaling limit. 



\F{L)) 

While Zi = {F[s) / F{y)Y^^^ violates the Slavnov-Taylor identity which instead 
would demand Zi = {F {s) / F [L))^^"^^ , it thus reproduces the correct scaling limit for 
the gluon propagator. This will be reflected in the numerical solutions, in particular, 
in the leading logarithmic behavior of the gluon propagator. 

Note that eqs. (83) and ( ^5| ) are equivalent to order in perturbation theory. 



One-loop scaling requires, however, that a certain class of diagrams of the pertur- 
bative series is subsummed, which includes contributions of orders g'^ and higher. 
It should therefore not be surprising that a truncation scheme neglecting 4-gluon 
correlations (which appear at order g'^) does not automatically reproduce the per- 
turbative anomalous dimensions. We have seen that this problem is also related to 
some violation of the implications of gauge invariance. 
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We will proceed mainly with discussing the solutions to (76) and (pBD analyti- 



cally as well as the numerical results, and use the alternative version (|8^) only for 
comparison and in order to demonstrate that our main conclusions such as the exis- 
tence of the infrared fixed point do not depend on the particular choice for Z\ in the 
3-gluon loop. We will verify that the sole effect of this choice is to fix the behavior 
of the solutions, in particular, the gluon propagator, at short distances where their 
anticipated asymptotic perturbative form was used as the guiding principle. 

7 Asymptotic Expansion and Scale Invariance 

The set of equations ( [76| ) and ( p5D does not depend on the scale s showing that 
its solutions represent gluon and ghost renormalization functions obeying one-loop 
scaling at all scales, i.e., = 0, see eqs. ( |69D and (^). We have thus rewritten the 
problem in terms of the renormalization group invariant functions F{x) and R{x). In 
particular, the scaling behavior of the propagators follows trivially from the solution 
for the non-perturbative running coupling F/Pq. 

Similar to our previous solution to the gluon DSE in Mandelstam approximation, 
a detailed analyses of the solutions in the infrared in terms of asymptotic series is 
necessary in order to obtain numerically stable iterative solutions. Due to the nature 
of the coupled set of equations a recursive calculation of the respective coefficients of 
the asymptotic series for F and R is considerably more difficult than in Mandelstam 
approximation in which simple recursion relations allowed to calculate these coeffi- 
cients to any desired order § . Fortunately, calculating the leading corrections to 
the asymptotic infrared behavior of R and F, as given above, proved sufficient to 
obtain numerically stable results. Thereby, for x < xq with some infrared matching 
point xq, the asymptotic series to at least second order, in a sense to be explained 
below, is used in obtaining iterative solutions for x > xq. The matching point xq 
has to be sufficiently small for the asymptotic series to provide the desired accuracy. 
However, limited by numerical stability, it cannot be chosen arbitrarily small either. 
This leads to a certain range of values of xq for which stable solutions are obtained 
with no matching point dependence to fixed accuracy. We verified that the additional 
inclusion of third order contributions in the asymptotic series has no effect other than 
increasing the allowed range for the matching point. 

The discussion of the solutions in the infrared is alleviated by the observation 
that one of the equations, eq. (^), 

R{x) _ ,f f 'dy p. X 1 



F^{3 



[ ^ R{y)F'''iy) - ^ R{x)F'-'ix)) , (88) 
Jo y ^ / 



can be converted in a first order homogeneous linear differential equation for R{x). 
Differentiating eq. ( p8[ ) with respect to x one obtains, 

R'{x) = T{x)R{x) , T{x) := + ^ " ^^') • (89) 
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The second equation to solve, eq. (|8l 
7x2 




3 R\x) 2Jo y \ R\y) 
where we have used that 

xl r%R(^y)F^ii^y) = -xl r^(RF^'-ba^V) + -^^^ + Ax 

8 A y2 yyj g y2 \ ^ y g i _ ^ 

with A = 1 ^(^RF^^ -ba'^^y^) . (91) 

From the leading infrared behavior, i.e., F — > a and R hx'^ for x — > 0, we see 
from eq. ( [8^ ) that an asymptotic infrared expansion of \/{R^F^^'^^) has to contain 
powers of x'^ as well as powers of x in subsequent subleading terms. This motivates 
the following Ansatz, 

S=Af 

R{x) = bx" Ci„„ x"^^+=^"'^+'(^+2«) (92) 

l,m,n=0 
T.=N 

F{x) = a J2 ^i^n a;'"^+3'^'^+'(i+2«) , with S := / + m + n (93) 

l,m,n=0 

and Cooo = -Dqoo = 1- The additional fractional power x^ in these expansions was 
introduced to find the most important subleading behavior possible from the consis- 
tency in the infrared. We will determine it to be ~ 2.05. With 2<3k<H-2k!^3 
this shows that for not too large orders N powers of different orders in this expan- 
sions do not mix in their successive importance at small x. Furthermore, we have 
repeatedly subtracted leading infrared contributions explicitly from integrals such as 
the one in ( |8^ above, assuming that the remaining contributions are integrable for 
X — > 0. For subleading contributions to R and F suppressed by powers of x'^ with 
u ~ 2.05 this is justified a posteriori. 

With the series ( p3[ ) we can calculate T(x) in eq. ( ^8|) to the same order N , 

r(x) = - ^^imn a;""""" , with Timn ■= mu + 3nK + 1{1 + 2k) , (94) 

i,m,n=0 

where the coefficients Eimn can be straightforwardly calculated from the coefficients 



Dimn of F (with -Eqoo = !)• The solution of eq. (88) for R with the integration 
constant set to b is then given by 



C TV 

R{x) = bx'^ exp < K 

I E=l 



Tlr. 
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Expanding this series to an appropriate order allows to relate the coefficients Cimn 
to Eijnn and thus to Dimn- For = 1 the result is, 

= rT^^ioo = (2a + 2K)+'^) ^100 

Coio = -^010 = ( - - ^ - ^ + ) ^010 (96) 

Com = g-E-ooi = (^3 ~ ~'~ ^ooi • 

At higher orders in this procedure recursively yields relations that uniquely de- 
termine the coefficients C in terms of the coefficients D. Analogous relations are 
obtained from eq. ( |90[ ) by expanding all ratios of R and F which occur with depen- 
dence on X and y, and by comparison of the respective orders, C'(x'^''""~^'^), on both 
sides. We verify from eq. ( |90| ) that to leading order, 0{x~'^'^): 

11 ^3 I 1 1 \ a'^^ 

T^' (97) 



62ai-2'5 \2 2-K 3 4^7 62 

which was used to determine k with a = /3o5c = ((9/^4) (1/k — l/2))~^ (see eqs. 
and (^^). At order = 1 we obtain. 



Oix"-^^) ■ — (Doio + 2(Coio - ^Z^oio)) = 

a 

0(x') : — (D„oi + 2(0„„i - SDm)) = -("'/(s) + 

a 

l(2T^+2^)-54)(^"'-*"»»'' 

0(1): — (Dioo + 2(Ci„o-*Oi„o)) = — + 

a a^" 

with /(k) := 7/(2(3 k)) - 17/(2(2 + k)) - 9/(8(1 + k)) + 7/k + 7/(8(1 - k)). 

These equations together with eqs. (|96| ) determine the coefficients D and C to 
lowest non-trivial order. In particular, we have 3 decoupled sets of 2 equations each 
for 2 of the constants, respectively. For {l,m,n) = (1,0,0), c./., eq. (|100|) , we obtain, 

Cioo ^ 0.05554 6^ A , and Dioo =^ -0.6992 6^ A . (101) 



The set of equations for {l,m,n) = (0,1,0) is homogeneous, see eq. (98). The deter- 
minant of its 2-dimensional coefficient matrix is zero for 



-6 - K - 3 k2 ± \ -4 (-26 - 23 k) k2 (3 + 2 k) + (6 + K + 3 ^2)2 

WT^) ■ 
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with one positive root for the plus sign which determines the positive exponent v. 
The coefficients Cqio, -Dqio fohow up to a common factor. With k = (61 - Vl897)/19 
we obtain, 

V ~ 2.051 and Cqio = -0.0124L»oio • (103) 

For (/, m, n) = (0, 0, 1) the scale of the coefficients is set by the inhomogeneity, b^f{K), 
in eq. (^) and we obtain, 

Cooi ^ 1.969 6^ , and L»ooi - -26.52 6^ . (104) 

Using these results, higher orders, though increasingly tedious, can be obtained re- 
cursively by analogous sets of equations. The general pattern is such that the order 
= 1 above fixes the scales for higher order coefficients. This allows us to define 
scale independent coefficients C and D by extracting their respective scales according 
to the exponent Timn = in^^ + 3n/t + /(I + 2k) of x for a given set (/, m, n), 

Clmn =■■ Cimny'^'^^h^ A' , and Amn =: A„„ 63^+24- , (105) 

where the scale of the powers of x'^ is set by t which we fix for convenience to 

t := -Dno (106) 

i.e., Cqio — 0.0124 and -Dqio = —1- We summarize the values of the coefficients C 
and I) for iV = 2 in table |. 



(/, m, n) 


(2,0,0) 


(1,1,0) 


(1,0,1) 


(0,2,0) 


(0,1,1) 


(0,0,2) 


C 


-0.1042 


-0.3034 


-7.933 


-0.2160 


-11.55 


-151.0 


D 


0.5246 


1.590 


40.10 


1.226 


60.98 


766.8 



Table 1: Coefficients of the asymptotic expansion for N = 2. 



The constant A as given in (^) is determined numerically in the iterative process. 
Of the remaining two parameters b and t in the asymptotic forms, one can be related 
to the overall momentum scale so that we are left with one independent parameter. 
Note that we neither specified the momentum scale a (in x = k"^ /(t) nor the infrared 
constant h so far. The problem is scale invariant, however, i.e., a change in the scale 
a according to a —> cr' = cj/A or, equivalently, x — > x' = Ax can be compensated by 

h ^ y = b/X'' , and t ^ t' = t/X" . (107) 

We can thus choose the scale without loss of generality such that the positive number 
6=1. The parameter t can in principle be any real number including zero. We can 
find numerically stable iterative solutions for not too large absolute values of t (see 
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Figure 6: The numerical solutions of F{x) and K{x) for t = and 6=1 together 
with their asymptotic expansions to order = 1 as well as = 2 at small x. 

below). Furthermore, it can be verified numerically, that a solution for a value of 

5 7^ 1 for fixed t is identical to a solution for 5 = 1 and t' = t}fl'^^ if x is substituted 
by x' = x/b^^'^. This is the numerical manifestation of the scale invariance mentioned 
above (for A = 1/6^/'^). Note that under scale transformations (|107D the constant 
A trivially transforms according to its dimension, A — > A' = A/X, without any 
adjustments from the way it is calculated, since 

A' = A/X = lunU r %R{y')F^\y') - b' a^' . (108) 

In figure ^ we plotted the numerically obtained solutions for F{x) and R{x) for 

6 = 1 and t = at small x together with their respective asymptotic forms to order 

= 1 and N = 2. The contributions of the order = 2 in the asymptotic expansion 
become comparable in size to the lower order at about x ~ 0.2 and dominant for 
increasingly higher values of x. The error in the asymptotic series being of the order 
of the first terms neglected, this is the usual indication for the range of x in which the 
asymptotic expansion to the given order can yield reliable results. In the particular 
calculation we used a value of about xq = 0.008 for the matching point of the iterative 
process to the asymptotic result. This is obviously well below the estimated range of 
the validity of the asymptotic expansion. 
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0.1 1 10 100 10" 10^ 10^ 

Figure 7: The numerical solutions of F{x) and R{x) with 6 = 1 for different values of 
the parameter t = {—4, —2, —1, 0, 1, 2, 4, 8, 16} (solid lines represent t = solutions). 

8 Numerical Results and Perturbative Limit 

The numerical results were obtained with the order = 1 in the asymptotic expan- 
sion for which no dependence on xq was observed for xq < 0.1. Higher values for the 
matching point are possible at higher orders. This shows that the expansion to next 
to leading order in the infrared {N = 1) is sufficient to find numerically stable results 
for a satisfactory regime of matching points. Similar results follow for t ^ 0. We cal- 
culated F and R for several values in the range —5 < t < 16. At lower negative values 
the procedure became numerically unstable due to a developing (tachyonic) pole in 
F{x). The fact that the integral equations for R and F possess a one-parameter 
family of solutions characterized by t is in fact the reason for the necessity of the 
infrared expansion up to next to leading order, since no stable solution can be found 
numerically without fixing the leading x-dependence of F{x) at small x by fixing 
the parameter t. This is a boundary condition to be imposed on the solutions from 
physical arguments. 

In figure |^ the numerical results are plotted for different values of the parameter 
t (all with 6=1). Perturbatively, we expect R{x) to approach a constant value and 
F{x) l/ln(Ax) for x — > oo. The reason we introduced the constant A in the one- 
loop running coupling F is that we fixed the momentum scale in our calculations 
by arbitrarily setting 6=1. The relation between the scale of perturbative QCD 
Aqcd and a cannot be determined this way. Therefore, we set Aq^j, = a/X for 
some scale parameter A. Qualitatively, all solutions display a similar behavior at 
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high momenta. The solutions for positive values of t seem to have more residual 
momentum dependence in R at high momenta than those for t < 0. A very good 
fit to the perturbative form can be obtained for t = —4.2. For negative values of t, 
however, the running coupling, as(/i) = F {s) / {An Po) , has a maximum, a,„ax > etc, 
at a finite value of the renormalization scale s = iJ-'^/cr. This is because the dominant 
subleading term of the running coupling in the infrared is determined by t, 

F{x) ail-tx" + Dooi x^'" + Dwo x'^^^") , x^O . (109) 

With ~ 2.05 < 3k < 1 + 2fi: and -Dqoi < 0, it is clear that for t < the running cou- 
pling increases for smallest scales close to /i = before higher order terms dominate. 
There necessarily has to be a maximum a,„ax > etc at some finite scale fi for any so- 
lution with t < 0. Such a maximum at finite scale implies a zero of the beta function 
which cannot be interpreted as a fixed point (since it appears at a finite scale fio). It 
is easy to see that the slope of the beta function P{g) at such a point ^^ax is infinite. 
This precludes its possible classification as being infrared or ultraviolet stable. The 
present case (t < 0) corresponds to a double valued P{g) with a positive branch and 
a negative branch coexisting for gc < 9 < 9max- In lack of an understanding of this 
scenario, we therefore disregard solutions for t < as likely to be unphysical. 

In contrast, for t > 0, Qc = a{fi = 0) is the only maximum of the running coupling 
for all real values of the renormalization scale, and ac is thus a true infrared stable 
fixed point. Comparing the behavior of the resulting gluon and ghost renormalization 
functions in the ultraviolet we observe that, for the t > solutions, the case t = 
yields the best resemblance of their one-loop anomalous dimensions. We therefore 
restrict the discussion of our results to the case t = from now on, interpreting the 
existence of solutions for t 7^ as a mathematical peculiarity. 

In order to obtain an estimate of the factor A to relate the infrared scale a to 
the perturbative QCD scale Aqcd, we compare the results for the renormalization 
functions Z{x) and G{x) (with x = k'^ /cr) to their respective one-loop forms. 



Vln(A;VA2^j; VMAx 



^ / ln(^VA^,,) \ ^ / ln(A.) 

[HkyAl^^)) lln(Ax)' • ^ ' 



In figure |8| we plot the functions Z(x) and G{x) as obtained from F(x) and R{x) 
(for t = 0) according to the parameterization ( |70|) along with their respective one- 
loop forms using A = 1.8. The renormalization point s = jj^ /c = 1.18 • 10'' (with 
Z[s) = G{s) = 1) was chosen to be the numerical ultraviolet cutoff of the particular 
calculation for aesthetic reasons. A different choice would result in no more than a 
constant vertical shift of the curves on the logarithmic plot, and the choice of the 
renormalization point is thus absolutely inessential for the present consideration. It is 
possible to match one of the renormalization functions Z and G with its perturbative 
form even closer than in this figure by varying the value of A in the perturbative 
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Figure 8: The renormalization functions Z{x) and G{x) obtained for t = and their 
leading logarithmic forms (with 7 —> 1 — 2(5 for the gluon, and j ^ 5 for the ghost). 

logarithms. This will, however, at the same time affect the perturbative form of the 
other renormalization function in the opposite direction, e.g., for A approximately 
4 ~ 5 the perturbative logarithm of the gluon is hardly distinguishable from Z(x) for 
a; > 10 on the scales as used in figure |^ while a similarly good fit of the one-loop 
form to G{x) is obtained for 0.2 ~ 0.3. To fit the perturbative logarithms at high 
momenta as close as possible to both renormalization functions with a unique value 
of A, we can restrict its value approximately to the range 1.5 < A < 2. 

For a scale factor in the same range, 1.5 < A < 2, the resemblance of the leading 
logarithms at high momenta can be optimized by an unphysical choice of the param- 
eter t of approximately t = —4.2. The according solutions approach their leading 
logarithmic behavior at momenta as low as k'^ ~ lOcr and are practically identical to 
the latter for all momenta higher than that. In contrast, as mentioned above, increas- 
ing positive values of t lead to solutions fit by the leading perturbative logarithms 
decreasingly well (see figure |9|) . 

To summarize the discussion of our results at high momenta, we conclude that we 
obtain quite good agreement with the leading logarithmic behavior (as known from 
perturbation theory) for the renormalization functions Z and G at high momenta. 
The perturbative QCD scale parameter Aqcd is thereby related to the infrared scale 
a by Aqqq = cr/A. The factor A can be numerically estimated to be approximately 
1.5 — 2. 

Note that even though we introduced a numerical ultraviolet cutoff, for sufficiently 
high values, our solutions to the ultraviolet finite equations for F{x) and R{x) show 
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Figure 9: The renormalization functions Z{x) and G{x) for t = {—4.2,0,16}. Their 
leading logarithmic form at large x is not given (see fig. |8|), it is indistinguishable 
from the t = —4.2 results for x > 10 (and with A in 1.5 — 2). 



no dependence on this cutoff. There are no ultraviolet boundary conditions imposed 
on the solutions. Fixing the infrared boundary condition by choosing a value for the 
parameter t (and the scale by setting b = 1) we have no further influence on the high 
momentum behavior of the solutions which selfconsistently results from the iterative 
process alone. Deviations from leading logarithmic behavior can have several origins: 

• Even though the gluon and ghost renormalization functions obey one-loop scal- 
ing, the non-perturbative running coupling, asifJ-) = F(s)/(47r/3o), contains 
two- and higher loop contributions (if certainly not all of them). Therefore, 
the fact that the results for t = approach the leading logarithmic behavior 
at higher momenta than some of the apparently unphysical results for negative 
values of t is not unreasonable. 

• Depending on t also, the function R{x) can have residual momentum depen- 
dence even at very high momenta. In figure the solution for, e.g., t = —4 
approaches the value one at momenta as low as k"^ ~ lOo", whereas for positive t 
considerable logarithmic momentum dependence can be left even at ~ 10^ a 
(see fig |7| for t = 16). While this has of course no influence on the running cou- 
pling, it does affect the asymptotic behavior of the renormalization functions. 

• Since ln(A2;) = InA -|- Inx, all different asymptotic forms used above yield an 
equivalent leading logarithmic contribution (in the limit x — > oo). The relevant 
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question is the one for the momentum scale above which the results should 
follow their leading logarithms in phenomenological applications. 

Fixing the momentum scale in our calculations from the rough estimate a = 
AAq(2J3 given above is not satisfactory. Not only that the parameter A is poorly 
determined, also the value of Aqqd, appropriate for the present subtraction scheme 
(and Nf = 0), is not a phenomenologically well known quantity. We discuss in more 
detail in the next section, how to avoid this problem by fixing a directly, independent 
of the precise values for A and thus Aqcd- As a result of this, we find that /c^ ~ lOo" 
corresponds to a momentum scale of IGeV, and we conclude that the high momentum 
behavior of the solutions for t = is in reasonable agreement with the perturbative 
result. 

In connection with the renormalization procedure discussed in the previous section 
we motivated the substitution of the cutoff dependence in the constant Zi by a 
dependence on the integration momentum to improve the ultraviolet behavior of the 
solutions. In figure ^ we compare the t = results for F{x) and R{x) obtained 
this way (from eq. ( |85|) for the running coupling F) with those of an analogous 
calculation using (|8^ ) with Zi = czj^ = 1 instead. It can be seen clearly that both 
solutions represent identical results in the infrared. This shows that the existence of 
the solutions and, in particular, of the infrared fixed point in the present truncation 
scheme is no artifact of the treatment we suggested for the renormalization constant 
of the 3-gluon loop Zi. Furthermore, it verifies that this treatment does indeed 
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Figure 11: The running coupling asip) for t = (left), and the corresponding /3- 
function (right) in comparison with its leading perturbative form (one-loop) and the 
/3-function as a result of the Mandelstam approximation from ref. Q. 

improve the ultraviolet behavior, noting that R{x) as obtained for Zi = 1 does not 
approach a momentum independent constant in the ultraviolet on any reasonable 
scale. 

9 Discussion of the Running Coupling 

The running coupling of the present scheme is given in eq. (^) by 

where s = fjp/cr involves the scale a related to the infrared boundary condition 6 = 1. 



The result for 05 (/i) is shown in figure 11. To assign a value to the scale a we use 



the phenomenological value of the running coupling as{Mz) = 0.118 at the mass 
of the Z-boson , Mz = 91.2 GeV For the t = solution we obtain this way 

Ml /a ~ 70000 and therefore, a ~ (350MeV)2. From the comparison of the solutions 
for the renormalization functions Z and G to their leading logarithmic form at high 
momenta, this corresponds to a perturbative scale Aqcd in 250 ~ 300 MeV ( for 
A in 1.5 ^ 2). This is not a very significant result, however, since the value for 
Aqcd is dependent on the number of quarks Nf and the order of the perturbative 
expansion. In particular, in the present framework we expect it to change when 
quarks are included. The determination of the momentum scale in our scale invariant 
calculation from the dimensionless function F{s) (by as{Mz) = 0.118) is independent 
of the value for the scale parameter A and thus of the value for perturbative scale 
Aqcd = \/cr/A. Equally independent of this value is the ratio of the Z~ to the r- 
mass, Mz/Mj- ~ 51.5, and we obtain from this ratio the value as{MT^) = 0.38 for 
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the running coupling at the mass of the r-lepton. Even though this is in compehing 
agreement with the experimental value [Q, we deliberately do not want to claim it to 
be a theoretical prediction. We regard this result as nothing but a consistency check 
of the present calculation which happens to yield better numbers than the present 
level of truncations and approximations might suggest. 

The right panel of figure 11 shows the corresponding /3-function with its two 
fixed points in the ultraviolet at g = and in the infrared at g = Qc ~ 10.9. For 
comparison we plotted the leading perturbative form for g ^ (and Nj = 0), 
P{g) — > —f3og^ = — ll/(167r^) g^, and the infrared singular non-perturbative result 
obtained from the analogous subtraction scheme in the Mandelstam approximation 
according to ref. |^] . The solution to the coupled system of gluon and ghost Dyson- 
Schwinger equations yields better agreement with the leading perturbative form at 
small g, since due to neglecting ghosts, the leading coefficient was obtained to be 
Po = 14/(167r^) in Mandelstam approximation j^. While this can be regarded as 
small quantitative discrepancy, the significant difference between the present result 
and the Mandelstam approximation for g > 7 once more demonstrates the importance 
of ghosts in Landau gauge, in particular, in the infrared. 



10 A Comparison to Lattice Results 

It is interesting to compare our solutions to lattice results available for the gluon 



propagator |34, ^] and for the ghost propagator |51] using lattice versions to imple- 
ment the Landau gauge condition supplemented by different procedures to eliminate 
Gribov copies (at least approximately). Recently, for the pure SU{2) lattice gauge 
theory, the influence of such copies of gauge equivalent conflgurations present in 
the conventional Landau gauge, has been systematically investigated for gluons and 
ghosts in 



In figure 12 we compare our solution for the gluon propagator to the data from 
ref. [35|. The momentum scale in our results, fixed from as{Mz), is not used as a free 
parameter. In order to account for the units used in ref. |^5|, we plotted the gluon 
propagator, normalized according to Z{x = 1) ~ 11.3, as a function of the momentum 
X = k'^a? in units of the inverse lattice spacing, using = 2 GeV corresponding to 
the value /3 = 6.0 for SU{?>) used in [^. According to the authors of |3^, the arrow 
indicates a bound below which finite size effects become considerable. As can be seen 
in figure |l^, the essential features of our solution are beyond the scope of present 
lattice calculations. Due to different normalizations, lattice sizes and values of the 
lattice coupling /?, it is not quite obvious that presently available lattice data for the 



gluon propagator in Landau gauge from the different groups p3] and |35] is indeed 



consistent. Since we prefer to adjust the units of our calculations rather than those 
of the lattice data, we apologize for choosing a particular set of data (from figure 3 
in ref. [^5|) and refer the reader to the literature, in order to asses the consistency of 
the different lattice results. 
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Figure 12: The numerical result for the gluon propagator from Dyson-Schwinger 



equations (solid line) compared to lattice data from fig. 3 in [35 



In figure 13 we compared our infrared enhanced ghost propagator with normaliza- 



tion such that G{x = 1) = 1 to the results of ref. [51|. We choose to display the lattice 



results for the symmetric 24 and the non-symmetric 16 x 32 lattices for SU (3) at 



P = 6.0 from figure 1 in ref. [51| up to x ~ 3. Identical results modulo finite size 
effects were obtained for an 8^ lattice (see ref. |5l|). Again, x = k'^o? with = 2 
GeV and the momentum scale in our results, fixed from the Z-mass as described in 
the last section, is not adjusted. 

The data extracted from the long direction of the 16^ x 32 lattice might indicate 
the existence of a finite maximum in the ghost propagator at very low momenta. 
The fact that the two lowest data points in this case lie significantly below their 
neighbors of the 24^ lattice was interpreted by the authors of ref. |51| as a genuine 



signal rather than a finite size effect. The reason for this being that, on the S'^ lattice, 
an enhancement due to finite size was observed for the lowest point in contrast to 
the shift downwards of the two points from the 32 direction of the non-symmetric 
lattice. No such maximum was observed on any of the smaller lattices. Our results 
do not confirm the existence of such a maximum in the ghost propagator but coincide 
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Figure 13: The numerical result for the ghost propagator from Dyson-Schwinger 
equations (solid line) compared to data from fig. 1 in ref. [^] for 24^ and 16^ x 32 
lattices, and a fit of the form c/x + d/x^ (with c = 0.744, d = 0.256) as obtained in 
ref. [^] for x >2 (dashed line). 



nicely with all those data points of the differently sized lattices that lie on a universal 
curve. In addition, the 24^ and 16^ x 32 lattices are of roughly the same size as the 
unsymmetric 16^ x 48 and 24^ x 48 lattices used for the gluon propagator in ref. 
||35|| . Considering their investigation of finite size effects due to deviations between 
different components of the gluon propagator at small momenta, one might be led to 
question the significance of the maximum in the ghost propagator observed for one 
particular set of data in ref. |51] at momenta too low to yield finite size independent 
results for the gluon propagator, c./., ref. [35| on even larger lattices.^ 

It is quite amazing to observe that our numerical solution fits the lattice data 
at low momenta [x < 1) significantly better than the fit to an infrared singular 
form with integer exponents, Dcik"^) = c/k"^ + d/k'^, as given in ref. |51]. Clearly, 
low momenta (x < 2) were not included in this fit, but the authors conclude that 
no reasonable fit of such a form is possible if the lower momentum data is to be 

Note the different a;-ranges used in figures |l| and |l| and the position of the arrow indicating 
finite size effects to become considerable. 
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included. Therefore, apart from the question about a possible maximum at the very 
lowest momenta, the lattice calculation seems to confirm the existence of an infrared 
enhanced ghost propagator of the form Dq ~ l/{k'^)^^'^ with non-integer exponent 
< K < 1. The same qualitative conclusion has in fact been obtained in a more 
recent lattice calculation of the ghost propagator in SU{2) |5^, where its infrared 
dominant part was fitted best by Dq ~ l/{k'^y^'^ for an exponent k of roughly 0.3 
(for (3 = 2.7). This also is in qualitative agreement with the SU{2) calculations of 



ref. |51], again, with the exception of one data point for the lowest possible lattice 
momentum. 

Furthermore, in refs. |51, 52] the Landau gauge condition was supplemented by 



algorithms to select gauge field configurations from the fundamental modular region 
which is to eliminate systematic errors that might occur due to the presence of Gribov 
copies.0 Thus, the good agreement of our result with the lattice calculations suggests 
that the existence of such copies of gauge configurations might have little efi^ect on 
the solutions to Landau gauge Dyson-Schwinger equations. This could also explain 
the similarity of our solutions to the qualitative behavior obtained by Zwanziger for 



gluon and ghost propagators from implications of complete gauge fixings |3£]. 

11 Confined Gluons 

Since we calculated Euclidean Green's functions, or Schwinger functions, for gluon 
and ghost correlations, we resort to the Osterwalder-Schrader reconstruction theo- 
rem which states that a Garding-Wightman relativistic quantum field theory can be 
constructed from a set of Schwinger functions if those Euclidean correlation func- 
tions obey certain conditions, the Osterwalder-Schrader axioms [^]. In particular, 
the axiom of reflection positivity for Euclidean Green's functions is a direct result 
of the positive definiteness of the norm in the Hilbert space of the corresponding 
Garding-Wightman quantum field theory. This condition involves arbitrary partial 
sums of n-point functions and is hardly provable in its general form. The special 
case of a single propagator G{x — y) reads, i.e., the lowest partial sum, 

j d^xd\j f{-xo,^)G{x-y)f{yo,y) > (113) 

where / G 5+(IR^) is a complex valued test (Schwartz) function with support in 
{(xojx) : xq > 0}. This special case of refiection positivity can be shown to be a 
necessary and sufficient condition for the existence of a Lehmann representation of the 
Wightman distribution (the Minkowski space propagator) corresponding to G{x — y) 



1 54]. Therefore, the construction of a counter example to this condition (by a suitable 
choice of /) is sufficient to prove that the Euclidean correlation function cannot 
represent physical particle correlations which can be interpreted as a manifestation 
of confinement. 



An investigation of the effectiveness of the stochastic overrelaxation method to achieve this is 
given in ref. [M 
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Figure 14: The one-dimensional Fourier transform of the gluon propagator, G{t, 



Heuristically, we can rewrite the condition (113) after three dimensional Fourier 
transformation as 



dt'dtfit')Gi-{t' + t),p)fit) > 0, 



(114) 



where G(xo,p) := / d^x G{xq,x.) e^P-^, and a separated momentum dependence of 
the analogous Fourier transform of the test function / has to provide for a suitable 
smearing around the three-momentum p. 

In figure |l^ we plotted this Fourier transform of essentially the trace of the gluon 
propagator, 



dpo Z{pI + p2 



27r 



+ p2 



(115) 



for p^ = and p^ = a /A. Choosing a sufficiently small p and a real function f{t) 
which peaks strongly around a point to/2 for which G(— to,p^) < 0, it is now very 
easy to construct a counter example to condition ( |113| ) for the gluon propagator. 
We interpret this result as confined gluon correlations described by the Euclidean 
Green's function Dfj^y. The ghost propagator in Euclidean momentum space, being 
negative for positive renormalization functions G{k'^) > which is the case at tree- 
level already, trivially violates reflection positivity. This is of course due to the wrong 
spin-statistic of the ghost fields. 
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Note that the existence of zeros in the Fourier transform (115) for sufficiently 
smah implies that the gluon propagator in coordinate space cannot have a well 
defined inverse for all (t,x). While the corresponding singularities in G~^(t,x) will 
cause problems in a Hubbard-Stratonovich transformation to bosonize effective non- 
local quark theories as the Global Color Model, see ref. [^], the same comment 
we gave at the end of section 4 applies here again. The infrared structure of the 
vertex functions, renormalization group invariance and the existence of zeros in the 
gluon propagator alone, all lead to the same conclusion: The combination of ghosts 
and gluons, g'^G'^{k'^)Z{k'^) /k'^ = as{k'^)/k'^ is the physically important quantity 
that determines the interactions of quarks in Landau gauge (and the current-current 
interaction of the Global Color Model). We have verified explicitly that the Fourier 
transform of as{k'^)/k'^ is free of such zeros and positive. 



12 Summary and Conclusions 

We presented a selfconsistent solution for the propagators of gluons and ghosts within 
the Dyson-Schwinger equation approach to QCD in Landau gauge without quarks. 
A direct comparison of this solution to the gluon propagator from the Mandelstam 
approximation to its Dyson-Schwinger equation shows the tremendous importance 
of non-perturbative ghost contributions which are neglected in the Mandelstam ap- 
proximation. 

We introduced a systematic truncation scheme to close the infinite tower of 
Dyson-Schwinger equations at the level of 3-point vertex functions. This was pos- 
sible by constructing these vertex functions from their respective Slavnov-Taylor 
identities. Such constructions, when contributions from irreducible 4-point scatter- 
ing kernels are neglected in the Slavnov-Taylor identities, allow to express all 3~point 
vertex functions in terms of the scalar functions which parameterize the propagators 
of gluons, ghosts and quarks. This procedure is consistent with the truncation of 
the Dyson-Schwinger equations in which irreducible 4-point correlations are also 
neglected. 

In this truncation scheme, at present without quarks, it is found that the leading 
infrared behavior of the gluon propagator is dominated completely by the ghost loop 
in the gluon DSE clearly demonstrating the importance of ghosts. In particular, an 
analytic discussion of the solutions to the closed set of equations for gluons and ghosts 
reveals that the infrared behavior of their 2-point functions is determined by an irra- 
tional exponent: The gluon propagator vanishes in the infrared as D{k) ~ (A;^)^'^^^, 
and the ghost propagator is infrared enhanced (as compared to a free massless particle 
pole), Dcik) ~ (l//c2)''+i, where k = {61 - Vl897)/19 ^ 0.92. 

The non-perturbative definition of the running coupling 05 (/u^) of our renormal- 
ization scheme is given by the renormalization group invariant product of the ghost 
and gluon renormalization functions G and Z, i.e., Airasifi'^) = g'^ G'^ {/j,^) Z (fj?) . 
Their infrared behavior implies the existence of an infrared stable fixed point at a 
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critical coupling, Oc — 9.5. 

Both, the gluon and the ghost propagator compare favorably with recent results 
from lattice gauge simulations. Since these studies implemented various algorithms 
to avoid possible Gribov copies, in particular, those of the ghost propagator, the 
compelling agreement with our solution for the ghost propagator indicates that the 
Gribov problem might have little influence on the solutions to Dyson-Schwinger 
equations. 

We have furthermore demonstrated that the gluon propagator obtained here vi- 
olates Osterwalder-Schrader reflection positivity, and hence has no Lehmann repre- 
sentation, which we interpret as a signal for confinement. To study quarks and, most 
importantly, their confinement, the present framework will be extended to include 
the quark propagator selfconsistently in future. 

We have given a variety of arguments in section 2 as to why our present results 
should not be sensitive to undetermined transverse terms in the 3-gluon vertex. 
However, the necessity of an additional assumption on the renormalization constant 
of the 3-gluon vertex to restore the ultraviolet behavior of the solutions, besides being 
due to neglecting 4-gluon correlations which appear at order g'^, might indicate that 
the possible presence of some additional purely transverse terms in the vertex should 
be taken into account also, in future. 

The conclusions of the present work rely on the fact that there are no such unde- 
termined transverse terms in the ghost-gluon vertex of the present truncation scheme. 
Therefore, in contrast to the ultraviolet (perturbative) behavior, the infrared behav- 
ior of the propagators, being entirely due to infrared dominant ghost contributions, 
is not affected by the details of the 3-gluon vertex. As for the gluon loop, we were 
able to proceed in close analogy to the Mandelstam approximation while accounting 
for ghost contributions. This allows a direct comparison and shows the influence of 
ghosts in Landau gauge clearly. 

Nevertheless, the construction of transverse terms by means of constraints arising 
from multiplicative renormalizability in the case of the vertices in QCD will be an 
important future project. Such a construction, available at present for massless 
quenched QED |29|, will be of particular importance when quarks are included in the 
selfconsistent framework. 

It will be furthermore important to compare our present results to calculations 
which do not rely on a one-dimensional approximation for the integrals. In or- 
der to overcome this approximation, it will be necessary to extend the technique of 
asymptotic infrared expansions to the original equations including their full angu- 
lar dependence. This analytic technique, developed and applied to the Mandelstam 
approximation previously [0, ^ and extended to the coupled system of gluon and 
ghost Dyson-Schwinger equations in the present paper, proved to be a necessary 
prerequisite to finding numerically stable solutions. 

Since confinement might be realized by a quite different mechanism in a different 
gauge, and in order to asses the possibility of gauge artifacts in our results, it will 
be important to study different gauges in parallel. Such a study of a coupled system 
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of equations, analogous to the present one, which determines the gluon propagator 
in the ghost-free axial gauge while accounting for its complete tensor structure, is 
currently in progress [p^. 
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Appendix 

A The Identity for the Ghost— Gluon Vertex 

Without irreducible ghost-ghost correlations, c./., eq. (p^, the Slavnov-Taylor iden- 
tity ( |T3|) decomposes into the irreducible ghost-gluon vertex function, G1^'^{x,y, z), 
and gluon and ghost propagators as follows. 



1 



A A dlD^^ix - u) Cf^iu, V, w) Dg{z - v)Dg{w - y) 



(116) 




(AidSj dSv d^D^^iy - u) GI^^{u, v, w) Dg{z - v)Dg{w - x) 



gr''DG{z-x)DG{z-y) . 



From the Slavnov-Taylor identity for the gluon propagator. 




and with the ghost propagator 



Dcix-y) = I A^DG(g)e^^(^-^) 



we obtain after Fourier transformation. 



ikpGf'^{k,q,p) 



Dg{p) DG{q) 



p2 



DG{-k)DG{q) 



gr'^ {2nf5\k + q - p) DG{-k) Dg{p). 



(117) 
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Here, G'^'^{k^q^p) denotes the ghost-gluon vertex in momentum space which corre- 
sponds to the conventions used in section 2 up to the total momentum conserving 
5-function, 



(fx<fy(fz Gf%x,y,z) e 



-ikx ^-iqy ^ipz 



{2T,f5\k + q-p)Gf'{p,q). 



Using this in eq. (|117D one obtains eq. (15) for the ghost-gluon vertex function. 



B 3— Gluon Vertex with Tree— Level Ghost CoupHng 

In this appendix we demonstrate that the Slavnov-Taylor identity for the 3-gluon ver- 
tex in presence of ghosts in Landau gauge has no solution assuming a tree-level form 
for the ghost-gluon coupling. Neglecting irreducible contributions from the ghost- 
gluon scattering kernel, i.e., setting G^^{p,q) = 5^^, the Slavnov-Taylor identity for 
the 3-gluon vertex function in Landau gauge simplifies to the following equation, 

ikpr^,p{p,q,k) = GiP) ^V^,{q)^^ - p^^(p)_^| . (118) 

Following the general Ansatz of ref. [^] for those parts of the 3-gluon vertex which 
can be constrained by its Slavnov-Taylor identity, we write, 

^fMupip, q, k) = -A{p^, q^; k^) i{p - q)p - B{p^, q^; k^) 6^„i{p + q)p 

-C{p^,q^;k'^){5^ypq - Puq,i) i{P - q)p + ^S{p^ ,q^ ,k^)i{ppqpku + p^qpkp) 

+ cyclic permutations . (119) 

The most general tensor structure contains these 10 terms plus 4 purely transverse 
terms involving 2 more scalar functions ||2^. The latter cannot be constrained by 
Slavnov-Taylor identities and are generally omitted. The scalar functions appearing 
in ( 119 ) have the following symmetry properties, 

A{x,y]z) = A{y,x-z) , B{x,y; z) = -B{y,x; z) , C{x,y; z) = G{y,x; z) , 

and S{x, y, z) is totally antisymmetric with respect to any two arguments. With this 
Ansatz in (118) we obtain the following equations from comparing the coefficients of 
the independent tensors on both sides: 

^i^u : (^A{x,y;z) - B{x,y;z)jy - (^A{x,y; z) + B{x,y; z)jx 

and -2B{x,y]z)+G{x,y]z){x-y) = , (121) 
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Q/iPu; 



-2A{z, y; x) + A{x, z; y) - B{x, z; y) 

--{z + x-y)S{x,y,z) = > 

same as above for q^q^ with x ^ y , 
-A{z, y; x) - B{z, y; x) + A{x, z; y) - B{x, z;y) =0 
-2B{x, y; z) - 2A(z, y; x) + 2A{x, z; y) - zS{x, y, z) 



The first observation here is that with (124) in ( |123| ) we find that 
S{x,y,z) = . 



(122) 

(123) 
. (124) 

(125) 



We could than in principle use ( |124| ) to eliminate B in (|122| ) thus obtaining an 
equation for A. Eq. ( 124| ) would then give us the corresponding function B which 
would fix C from ( |121 ). However, eliminating B from ( |122| ) yields, 



A{x, y; z) - A{x, z\ y) + A{z, y; x) 



Gjz) 

z{y) ■ 



(126) 



From the necessary symmetry of A{x, y; z) it follows that the l.h.s. in ( |126| ) is sym- 
metric under exchange of x and z whereas the r.h.s. is not. It is therefore not possible 
to find a solution A for this equation with the correct symmetry. 

This can be cured in two ways. The first is to allow for a non-trivial ghost-gluon 
scattering kernel of the form given in eq. (|lO| ) in section 2 which would not affect the 
ghost-gluon vertex since the correction to the kernel is transverse in the momentum 
of the incoming ghost. With the Ansatz (^) we can redo the above analysis obtaining 
modifications to the different terms on the r.h.s. of eqs. ( |12[1| ) to ( |122| ). In particular, 
eq. (|120|) is replaced by 



[A{x,y;z) - B{x,y;z)jy - (^A{x,y; z) + B{x,y; z)jx 

1 + ya{z,y;x 



\z{y) 

eq. (|123) by, 

- 2B{x, y; z) + C{x, y; z) {x - y) 



Z{x 

ya{z,y;x) xa{z,x;y) 



(127) 

(l + xa{z,x;y)j j G{z) , 



and eq. (|129D by 

- 2A{z, y; x) + A{x, z; y) - B{x, z; y) 



z{y) 

G{z) 



Z{x) 



G{z), (128) 



z{y) 



1 + ya(z,y;x)) . 



(129) 



Eqs. ( |123| ) and (|l2j) remain unchanged and, in particular, 5 = 0, which we used in 
(|129D already. Defining 



/(z, X] y) := ^p- (l + y a{z, y; zfj , 



z{y) 



(130) 
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we realize that the program outlined above to obtain the functions A, B and C can 
give a solution if f{x,y;z) = f{y,x]z). The equation to obtain A which replaces 
(126) is now, 

Mx, y, z) - A{x, z] y) + A{z, y; x) = f{z, x; y) . (131) 

Note that using ( |124| ) here again, this is yields, 

y, z) - B{x, y; z) = f{z, x; y) , (132) 

showing that a solution to ( |131| ) for A with B obtained from (|l2^) also satisfies 
(p^. With C from (|12|) we find that all necessary equations (|l2^,|l2|Jl2|) and 
(123, 124| ) are solved (with 5 = 0). So far, we did not specify the unknown function 



a{x,y;z) resulting from a non-trivial ghost-gluon scattering kernel. The necessary 
symmetry of f{x, y; z) provides certainly not enough information to fix a{x, y; z) 
completely. However, since we do not want to modify the scaling properties of the 
vertex artificially, it may seem reasonable to solve the system of equations with an 
Ansatz for f{x,y;z) as the most general sum of terms representing ratios G/Z with 
all possible arguments x, y and z obeying the necessary symmetry of / in its first two 
arguments. This restriction, which yields terms of the same structure as the r.h.s. of 
eq. ( |126| ) now in (|131| ), allows to determine the function a{x, y; z) up to an unknown 
constant A, 

._\(G{z^_G{y)\ _ A \G{y)(Z{y) , Z{y) ^\ 
"^"^'^'^^-y U(x) G{x)) 2y\G{x)\z{z)^ Z{x) 

The solutions for A, B and G follow as outlined above, in this case. 



2 \Z{x) Z{y) J 2 \Z{x) Z{y) 

1-A f G{x) _ G(y) \ _ A /G(x) _ G{y) \ 
2 \Z{x) Z{y)) 2 \Z{z) Z{z)) 



(135) 



x-y \Z{x) Z{y), 

This solution to the 3-gluon Slavnov-Taylor identity corresponds to a specific form 
for the ghost-gluon scattering kernel, in the simplest case, for A = 0, the dependence 
on ratios of gluon renormalization functions Z disappears and it is given by 

G^u{p,q) = V + ^ [Ofiiypq - Puqti) ■ (137) 
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The second way to account for the presence of ghosts in the Slavnov-Taylor identity 
in Landau gauge is to consider a possible dressing of the ghost-gluon vertex function 
as described in section 2. For the present truncation scheme this led to a vertex as 
given in eq. (|2^) from the new Slavnov-Taylor identity (|l3|). The structure of the 
corresponding ghost-gluon scattering kernel as far as it is relevant to the 3-gluon 
Slavnov-Taylor identity, however, is then simply 



Gf,u{p,q) = -^S^, , (138) 



which, used in the Slavnov-Taylor identity (P), led to eq. (|2^). The solution to this 
equation, as given in eqs. (25/2^) in section 2, can be obtained from an entirely 
analogous procedure as outlined in this appendix, setting, 

G{x)G{y) , , 



as it corresponds to the AbeJiaii-like identity (P4|). This results in the solution 
A{x,y\z) = A+{x,y;z), B{x,y;z) = A^{x,y;z) and C{x,y;z) = 2A_{x,y;z)/{x-y) 
with the notations of section 2. The latter way to implement ghost contributions in 
the construction of the 3-gluon vertex function from its Slavnov-Taylor identity is 
considerably simpler than the alternative as described in this appendix (based on a 
tree-level ghost-gluon vertex), and it is furthermore much more consistent with the 
assumption underlying the truncation scheme of section 2. 



C Functions in the 3— Gluon Loop 

Below we give the functions A^i and N2 which appear in the 3-gluon loop of the gluon 
DSE using the solution (25/^) for the 3-gluon vertex function and contracting the 
gluon DSE with the projector lZ^y{k) of eq. (p7[) (see sec. 4), 



Ni{x,y]z) 



29x x^ . 29y y 



Ay 



^ „^ ^ (9x2 + 50xy + 9y2) ^ 



Ax 



Axy 



9 {x + y) z' 
Axy 



Axy 



^2Ax^ -Wxy + y^ ^2Ay' 



(140) 
10 x y + 



2(z-x) 



and 



N2{x,y;z) 



5 x'^ + 41 x^ y + 5 X y2 



4x {y — x) 



+ 



x^ - lOxy + 24y2 
2 (2/ - ^) 



(141) 



^ x^ + 9 x^ y — 9 X y2 — ^ (2 x^ + 11 x y — 3 y^) z ^ (x + y) 



X z 



2x (x - y) 



4x (y — x) 
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